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EXECUTIVE  SUMMARY 


V  Fully  developed,  laminar  liquid-metal  flows,  currents,  and  power  losses  in  a  rectangular  channel  in  a 
uniform,  skewed  high  external  magnetic  field  were  studied  for  high  Hartmaim  numbers,  high  interaction 
numbers,  low  magnetic  Reynolds  numbers,  and  difleient  aspect  ratios.  The  channel  has  insulated  side  walls 
that  are  skewed  to  the  external  magnetic  field.  Both  the  perfectly  conducting  moving  top  wall  with  an  exter¬ 
nal  potential  and  the  stationary  perfectly  conducting  bottom  wall  at  zero  potential  act  as  electrodes  and  are 
also  skewed  to  the  external  magnetic  field.  A  solution  is  obtained  for  high  Hartmaim  numbers  by  dividing 
the  flow  into  three  core  regions,  connected  by  two  free-shear  regions,  and  Hartmaim  layers  along  all  the 
channel  walls.  Mathematical  solutions  are  presented  in  each  region  in  terms  of  singular  perturbation  expan¬ 
sions  in  negative  powers  of  the  Hartmaim  number. 

Numerical  calculations  are  presented  for  the  total  current  carried  by  the  core  region  between  top  and 
bottom  electrodes,  Joulean  and  viscous  power  losses,  and  chaimel  resistance  at  different  skewed  external 
magnetic  field  angles.  With  the  high  external  magnetic  field,  the  current  through  the  central  core  region  be¬ 
tween  the  electrodes  must  be  parallel  to  the  external  magnetic  field  lines.  The  two  side  core  regions  cany  no 
current  to  zeroth  order.  The  two  free-shear  layers  carry  less  current  than  the  central  core  region.  Both  free 
shear  layers  are  treated  rigorously  and  in  detail  with  fundamental  magnetohydrodynamic  theory.  The  solution 
for  the  free-shear  layer  velocity  profiles  are  solved  in  terms  of  a  complex  integral  equation.  The  integral 
equation  is  solved  by  the  method  of  quadratures  to  give  the  velocity  profiles,  viscous  dissipation  and  Joulean 
losses  in  the  free  shear  layers. 

In  addition  expressions  for  the  viscous  dissipation  in  the  six  Hartmann  layers  are  presented.:  The  best 
approximation  to  the  viscous  disssipation  in  the  channel  is  the  sum  of  the  0(M3/2)  contributions  from  the 
two  free  shear  layers,  the  0(M3/2)  contributions  from  the  two  Hartmann  layers  separating  the  free  shear 
layers  from  the  insulators,  and  the  0(M)  contributions  from  three  of  the  Hartmann  layers  separating  core 
regions  from  the  walls.  The  best  approximation  to  the  Joulean  power  losses  in  the  chaimel  is  the  sum  of  the 
0(M2)  contribution  from  the  centr^  core  region  which  carries  an  0(1)  current  between  the  electrodes  and  the 
0(M3/2)  contributions  from  the  free  shear  layers.  The  expressions  for  the  viscous  dissipation  and  Joulean 
losses  in  each  region  involved  the  products  of  universal  constants  electrical  potentials  and  geometric  factors. 
The  theoretical  magnetohydrodynamic  model  presented  here  was  developed  to  provide  data  to  help  in  the 
design  of  liquid-metal  current  collectors. 
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Fully  developed,  laminar  liquid-metal  flows,  currents,  and  power  losses  in  a  rectangular 
channel  in  a  uniform,  skewed  high  external  magnetic  field  were  studied  for  high  Hartmann 
numbers,  high  interaction  numbers,  low  magnetic  Reynolds  numbers,  and  different  aspect 
ratios.  The  channel  has  insulating  side  walls  that  are  skewed  to  the  externa!  magnetic  field. 

Both  the  perfectly  conducting  moving  top  wall  with  an  external  potential  and  the  stationary 
perfectly  conducting  bottom  wall  at  zero  potential  act  as  electrodes  and  are  also  skewed  to  the 
external  magnetic  field.  A  solution  is  obtained  for  high  Hartmann  numbers  by  dividing  the 
flow  into  three  core  regions,  connected  by  two  free-shear  regions,  and  Hartmann  layers  along 
all  the  channel  walls.  Mathematical  solutions  are  presented  in  each  region  in  terms  of  singular 
perturbation  expansions  in  negative  powers  of  the  Hartmann  number.  The  free-shear  layers  are 
treated  rigorously  and  in  detail  with  fundamental  magnetohydrodynamic  theory.  Numerical 
calculations  are  presented  for  the  total  current  carried  by  the  core  region  between  top  and 
bottom  electrodes,  Joulean  and  viscous  power  losses,  and  channel  resistance  at  different 
skewed  external  magnetic  field  angles.  With  the  high  external  magnetic  field,  the  current 
through  the  central  core  region  between  the  electrodes  must  be  parallel  to  the  external 
magnetic  field  lines.  The  two  side  core  regions  carry  no  current  to  the  zeroth  order.  The  two 
free-shear  layers  carry  less  current  than  the  central  core  region.  The  theoretical 
magnetohydrodynamic  model  derived  here  was  developed  to  provide  data  to  help  in  the  design 
of  liquid-metal  current  collectors. 


l^NTROOUCTiON 

'^Advanced  homopolar  electricaltMOTimery  is  being  con¬ 
sidered  for  many  new  applications^^his  applied  research 
requires  the  development  of  mechanisms  for  transporting 
high  current  at  low  voltages  between  rotating  and  stationary 
members  of  the  machines  with  minimal  losses  and  maximal 
operational  stability.  Thus,  liquid  metals  rather  than  more 
conventional  brush  technology  are  often  used  in  the  sliding 
electric  contact  region.  Generally,  these  machines  have  large 
external  magnetic  fields  in  the  region  of  the  current  collec¬ 
tor,  which  has  a  definite  effect  on  the  liquid  metal  by  Lorentz 
force  interactions.  Thus,  to  obtain  current  collector  design 
parameters,  one  is  interested  in  liquid-metal  channel  flow 
problems  with  applied  external  magnetic  fields  and  bound¬ 
ary  conditions  containing  combinations  of  moving  and 
fixed,  conducting,  or  insulating  walls.  These  boundary  con¬ 
ditions  then  roughly  correspond  to  a  rotor,  stator,  or 
sidewall  of  the  current  collector  and  the  liquid  metal  to  th^' 
sliding  electric  contact.  While  many  publications  trranhro- 
retical  magnetohydrodynamic  channel  duct  flow^^only  a 
very  small  percentage  of  these  investigation^ have  studied 
problems  with  moving  conducting  walls.^'Brown,  Reilly, 
and  Sondergaard"  have  studied  fully  develorod,  viscous  liq¬ 
uid-metal  flows  and  power  losses  in  a  rectangular  channel 
with  a  moving  conducting  wall  in  a  uniform, /xtemal  trans¬ 
verse  magnetic  field  for  moderate  Hartmanii  numbers.  The 
solutions  to  the  problem  were  presented  in  mfinite  series  of  a 
type  discussed  by  Hughes  and  Young.'/ The  rectangular 


channel  studied  had  insulating  side  walls  parallel  to  the  mag¬ 
netic  field  and  a  perfectly  conducting  moving  top  wall  and 
stationary  bottom  wall  perpendicular  to  the  field.  The  objec¬ 
tive  of  this  paper  is  to  treat  rigorously  and  comprehensively 
this  same  problem  with  two  major  changes.  The  homogen¬ 
eous  transverse  external  magnetic  field  is  changed  to  a  high 
fipld  strength,  homogeneous,  skewed  magnetic  field  and  an 
external  potential  is  applied  to  the  conducting  moving  wall, 
thus  making  the  channel  transport  current.  These  changes 
complicate  the  channel  problem  by  removing  symmetry. 
This  paper  should  be  considered  as  part  two  to  the  above- 
mentioned  paper. ' '  Alty presents  solutions  for  a  channel  j 
with  all  fixed  walls  and  a  skewed  magnetic  field.  His  core 
solutions  are  correct  but  his  free-shear  layer  solutions  are  i 
incorrect.  To  the  authors’  knowledge,  no  other  reference  in  j 
the  literature  addresses  this  particular  problem.  i 

The  purpose  of  this  paper  is  to  present  comprehensive 
two-dimensional  analytic  mathematical  solutions  for  the 
flow  velocity  and  electric  potential  in  accordance  with  the 
principles  of  singular  perturbation  th,ory  and  matched 
asymptotic  expansions  in  the  different  regions  of  the  chan¬ 
nel.  For  high  Hartmann  numbers,  the  channel  is  divided  into 
three  core  regions  ( a  central  core  and  two  side  core  regions ) , 
two  free-shear  layers,  and  Hartmann  layers  along  all  the 
walls.  These  solutions  enable  one  to  investigate  the  influence 
on  the  velocity  and  electrical  potential  of  insulating  side 
walls  and  of  an  external  skewed  high  magnetic  field.  The 
results  will  enable  the  derivation  of  the  equations  for  the 
quantities  of  technical  interest.  These  quantities  are  load 
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currents,  Joulean  and  viscous  power  dissipation,  and  chan¬ 
nel  electrical  resistance  as  functions  of  tne  skewed  external 
field  angle.  The  rectangular  channel  problem  approximately 
corresponds  to  an  asymmetric  liquid-metal  current  collector 
whose  channel  dimensions  are  small  compared  to  the  radius 
of  curvature  of  the  machine.  In  this  paper  a  skewed  magnetic 
field  was  chosen,  as  opposed  to  the  simpler  case  treated  ear¬ 
lier  of  a  transverse  magnetic  field,  because  many  current  col¬ 
lectors  have  skewed  magnetic  fields  in  relationship  to  the 
walls. 

In  this  work  analytical  expressions  are  derived  for  the 
load  current  between  the  electrodes  (i.e.,  moving  and  sta¬ 
tionary  conducting  wall),  Joulean  and  viscous  dissipation, 
and  channel  resistance  with  skewed  angle.  The  parameters 
of  variation  areas  follows;  applied  skewed  external  magnetic 
field,  fluid  electrical  conductivity  and  viscosity,  and  channel 
aspect  ratio.  In  the  work  presented  herein  the  pressure  gradi¬ 
ents  are  assumed  to  be  zero.  Numerical  results  in  nondimen- 
sional  form  are  presented,  along  with  the  general  analytic 
solutions  which  can  serve  as  benchmarkers  for  various  nu¬ 
merical  computation  procedures. 

This  paper  presents  the  first  correct  mathematical  treat¬ 
ment  of  the  free-shear  layers  resulting  from  skewed  external 
magnetic  field  interacting  with  rectangular  liquid-metal- 
filled  channel  described.  In  treating  the  magnetohydrodyna¬ 
mic  flow  due  to  electric  current  between  two  circular  elec¬ 
trodes  embedded  in  two  parallel  insulators,  with  a  transverse 
magnetic  field,  Hunt  and  co-workers  correctly  treat  a 
similar  free-shear  layer,  although  their  approach  is  quite  dif¬ 
ferent  from  the  present  one. 

A  number  of  steps  here  are  standard  parts  of  analytical 
solutions  for  fully  developed  magnetohydrodynamic  duct 
flows  at  high  Hartmann  numbers.''"'''  We  have  reproduced  a 
number  of  these  standard  steps  here  because  these  solutions 
are  not  widely  known  outside  the  relatively  small  magneto¬ 
hydrodynamic  community. 

II.  RECTANGULAR  CHANNEL  CONFIGURATION 
A.  Discussion  of  probiem 

The  problem  discussed  in  this  paper  consists  of  a  rectan¬ 
gular  channel  filled  with  a  liquid  metal  in  a  uniform,  skewed, 
high  external  magnetic  field'"  (see  Fig.  1).  The  Hartmann 
number  M  is  much  greater  than  one  1 ).  The  thin  per¬ 
fectly  conducting  top  wall  with  an  external  mean  potential  is 
assumed  to  move  at  a  velocity  component  C/,„  while  the  thin 
perfectly  conducting  bottom  wall  is  stationary  at  zero  poten¬ 
tial.  These  conducting  walls  are  sometimes  referred  to  as 
electrodes.  The  side  walls  of  the  channel  are  insulators.  We 
will  consider  the  case  of  fully  developed  laminar  flow  of  an 
incompressible  fluid  with  uniform  electrical  conductivity  a, 
and  uniform  viscosity//  ^ .  The  rectangular  duct  has  a  height 
of  L  and  a  width  oilaL.  We  also  assumed  that  no  secondary 
flows  are  generated  during  this  process,  that  there  is  no  vari¬ 
ation  of  the  duct  cross  section,  and  that  there  is  no  distortion 
of  the  external  magnetic  field  by  the  flow  (low  magnetic 
Reynolds  number). 

The  angle  between  the  skewed  external  field  B„  and  the/' 
axis  is  0.  In  Fig.  2  a  rotated  coordinate  system  (x,  77,  is 
shown  with  respect  to  the  coordinate  system  {x,y,z).6K  the 


FIG  I .  Rectangular  channel  with  skewed  external  magnetic  held  and  mov¬ 
ing  and  stationary  perfect  conductor. 


angle  between  they  axis  and  rj  axis  measured  positive  in  the 
clockwise  direction.  The  three  core  regions  Cl,  C2,  and  C  3 
are  also  shown  in  the  figure,  b  is  the  unit  vector  in  the  direc¬ 
tion  of  the  external  field. 

B.  Nondimensional  magnetohydrodynamic  equations 

The  nondimensional  magnetohydrodynamic  equations 
for  fully  developed  laminar  duct  flow  in  an  external  magnet¬ 
ic  field  at  an  angle  6  with  the  vertical  can  be  expressed  as 


W“'(vV)v  =  -  VP -f  jxb -f  Af  (la) 

j=  -V(i-t-vXb,  (lb) 

V-v  =  0,  ( Ic) 

V-j==0.  (Id) 


In  Eq.  ( 1  a ),  =  crP  0  Z,  //»  C/q  represents  the  interaction  par- 
ameter,the  ratio  of  pondermotive  force  to  the  inertial  force. 
M  =  LB„^(t/H f  is  the  Hartmann  number,"'  the  positive 
square  root  of  the  ratio  of  the  pondermotive  force'"  to  the 
viscous  force.  It  is  assumed  during  this  work  that  the  mag¬ 
netic  Reynolds  number'"  —  UJL  afi  (ratio  of  induced 
magnetic  field  to  external  magnetic  field)  is  <  1.  The  nondi¬ 
mensional  variables  in  Eqs.  ( la)-(  Id)  are  defined  as 


moving  perfect  conductor 


FIG  2  Rectangularchannelwithcoordinatesystems(x,j',z)and(x,77,|') 
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V  =  fluid  velocity  vector  normalized 

by£/o  =  V/C/o.  (2a) 

P  =  pressure  normalized 

by  orUoBlL  =  P*/aUoBlL,  (2b) 

j  =  electric  current  density  vector 

normalized  by  aU^  =  i/aUo  B„,  (2c) 

<f>  =  electric  potential  normalized 

by  (/o5oI  =  (^V(/o5„Z„  (2d) 

V  =  gradient  normalized  by  Z.  =  Z,V, 
b  =  unit  vector  parallel  to  the  external 


magnetic  field  =  cos  +  sin  ^  (Bo  =  fiob).  (2e) 
x,y,z  =  unit  vectors  in 

Cartesian  coordinates  (see  Fig.  1 ).  (2f) 

In  the  Cartesian  system  of  coordinates  ( see  Fig.  1 )  thej>, 
z  coordinates  are  normalized  by  L.  The  moving,  perfectly 
conducting  wall  with  velocity  U^^x  is  at  =  1;  the  stationary 
conducting  wall  is  at  >'  =  0.  The  insulating  wails  are  at 
z=  ±  a.  The  aspect  ratio  of  the  channel  is  represented  by 
la.  The  magnitude  of  the  external  magnetic  field  B„  and  the 
angle  6  between  the  field  and  the  vertical  is  represented  in 
terms  of  the  field  components  as 

B„=(Bl  +  B\)'‘\  (3a) 

0  =  arctan(5j/5^).  (3b) 

Equation  (la)  is  the  nondimensional  Navier-Stoke’s 
equation  with  an  external  jxb  magnetic  force.  Also,  Eq. 
( lb)  is  the  usual  nondimensional  expression  for  the  current 
density  in  the  channel  (ohms  law),  and  Eq.  ( Ic)  is  the  non- 
dimensional  expression  for  the  incompressibility  of  the  fluid. 
Equation  ( 1  d )  is  the  equation  for  the  conservation  of  charge. 

We  define  a  rotated  nondimensional  Cartesian  coordi¬ 
nate  system  (x,  r],  §)  with  the  rj  axis  parallel  to  the  external 
field  B„.  Therefore,  q  equals  the  unit  vector  b  (see  Fig.  2). 
The  coordinate  ^  is  shown  in  the  figure  where  |  is  the  unit 
vector  in  the  direction  Also  note  that  x  is  the  unit  vector  in 
the  axial  direction  along  the  channel  into  the  plane  of  the 
figure.  Therefore,  in  our  rotated,  right-handed  coordinate 
system  the  equations  describing  the  locus  of  the  top,  bottom. 


left,  and  right  walls  are,  respectively, 

77  =  sec  0 -t- (tan  5)1',  (4a) 

where  —acosd  —  sin  0<^<a  cos  6  —  sin  6, 

7]  —  tan  9^,  (4b) 

where  —  a  cos  cos  6, 

77=  —  a  csc9  —  (ctn  9)^,  (4c) 

where  —  a  cos  9  —  sin  0<|'<  —  a  cos  9, 

77  =  a  CSC  0  —  ctn  9^,  (4d) 


where  a  cos  9  —  sin  9<^^a  cos  9.  These  relationship.',  were 
derived  from  geometric  considerations. 

It  should  be  noted  that  the  unit  vectors  in  the  two  sys¬ 
tems  have  the  relationships: 

z  =  sin  01] -f  cos  0|,  (5a) 
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y  =  cos0i]  — sin0|.  (5b) 

In  our  rotated  coordinate  system  (x,  77,  ^)  the  physically 
realizable  variables  for  fully  developed  laminar  flow**  can  be 
expressed  as 

V  =  u(77,|)x,  (6a) 

P  =  P{r},^),  (no  axial  pressure  gradient),  (6b) 

4>  =  ^{r}4),  (6c) 

}  =  +  (6d) 


Substituting  these  variables  in  Eqs.  ( la)-(  Id)  results  in  the 
following  convenient  system  of  magnetohydrodynamic 
equations  in  terms  of  the  velocity  and  electrical  potential. 
The  electrical  potential  and  flow  velocity  are  variables  that 
are  measured  in  experimental  work  where 


d-u\ 

di  .  d<b  , 
A,  =  -  — .  A  =  -  37  + 

077 

^  +  4.  =  o. 

077  0|- 


(7a) 

(7b) 

(7c) 


We  must  now  consider  the  moving  perfectly  conducting 
top  wall  in  regard  to  the  current  density  j„,  and  potential  . 
The  nondimensional  current  density  in  the  wall  is 

j„,  =  -F  I),  (8a) 

where 

|  =  v„,xb  =  xXil.  (8b) 

Since  the  wall  is  a  perfect  conductor,the  conductivity  of  the 
wall  cTu,  ><7,  where  tr  is  the  conductivity  of  the  conducting 
fluid  in  the  channel.  For  ohms  law  in  the  magnetohydrody¬ 
namic  approximation,  it  can  be  shown  as  au,  —  00  that 
must  equal  ^  for  a  finite  current  in  a  perfect  conductor.  Max¬ 
well’s  equation  V  X  H  =  j  determines  the  induced  current 
in  the  system.  Therefore,  in  the  top  wall  +  where 

1^0  is  the  electrode  potential  at  ^  =  0  and  V  =  (0/ 
dr])i\  +  (0/0|^)|.  Similarly,  in  the  bottom  wall  =  0  and 
--  0  in  the  bottom  conducting  wall  of  the  channel. 

It  is  now  necessary  to  specify  the  boundary  conditions 
on  the  four  walls  of  the  rectangular  channel.  The  top  moving 
perfectly  conducting  wall  of  the  channel  has  the  following 
boundary  conditions  on  the  nondimensional  velocity 
and  the  potential  =  (j)-. 

u(.rj,^)  =  1,  (9a) 

=  <^0  + 1"  at  77  =  sec  0  -F  (tan  9)i,  (9b) 

where 


—  a  cos  0  —  sin  0<^<a  cos  0  —  sin  0. 

The  bottom  perfectly  conducting  wall  has  the  boundary  con¬ 
ditions 

u[7j,^)=0,  ^  =  Oat77=  (tan0)|,  (9c) 

where 


—  a  cos  0<^<a  cos  0. 

The  left-hand  insulating  wall  has  the  boundary  conditions 
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u{v4)  =  0, 

sin  6 y,  +  cos  0y^  =  0,  at  7/  =  —  a  esc  0  —  (ctn  6)^, 

(9d) 

where 

—  a  cos  0  —  sin  —  a  cos  6. 

Here  it  must  be  remembered  that  the  current  density  lines 
must  be  parallel  to  the  insulating  wall.  Similarly,  the  right- 
hand  insulating  wall  has  the  boundary  condition 

u(j7,^)=0,  (10a) 

sin  0 y,  -t-  cos  6 ji^  =  0,  at  rj  —  a  esc  6  —  ctn  6^  ( 10b) 

for  a  cos  d  —  sin  cos  0. 

In  this  work  the  Hartmann  number  M  is  always  as¬ 
sumed  to  be  high  (i.e.,  1 ).  In  this  case  it  can  be  shown 

that  distinct  subregions  of  flow  exist  in  the  channel  and  are 
shown  in  Fig.  3.  In  Appendix  A  is  a  derivation  of  the  order  of 
thickness  of  the  various  regions  in  the  channel.  The  angle  6 
has  the  range  of  values  0  <  0  <  arctan(2fl). 

In  the  inertialess  core  regions  Cl ,  C  2,  and  C  3,  the  first 
partial  derivatives  d  /d^  and  d  /drj  are  of  0(  1 ) .  In  the  inertia¬ 
less  Hartmann  layers  A  1-A  6,  the  layers  are  of  approximately 
0{M  ~ ' )  thickness  and  have  0{M)  normal  derivatives.  /, 
and  /j  are  inertialess  free-shear  layers  (or  interior  regions) 
which  have  the  following  characteristics: 

( 1 )  The  layers  separate  the  different  core  regions. 

(2)  The  layers  lie  along  magnetic  field  lines  through  the 
corners  and  are  driven  by  current  singularities  at  the 
corners. 

( 3 )  The  layers  have  thickness  ( see  Appendix 
A  for  detailed  derivation). 

(4)  The  first  partial  derivatives  are  of  orders  8/ 
(9J- =  0(  A/ '^-),  while  5/^77  =  C(  1 ). 

III.  ZEROTH  ORDER  MATHEMATICAL  SOLUTIONS  IN 
CORE  REGIONS 
A.  Development  of  theory 

The  generalized  mathematical  solutions  in  the  cores, 
where  the  magnetic  field  angle  6  has  the  limits  Q<d< arc- 


FIG  3.  Subregions  in  magnetohydrodynamic  channel  at  high  external 
magnetic  fields 


tan  (2a),  can  be  obtained  from  Eqs.  (7b)  and  (7c)  iftheterm 
ofO(A/"^)  is  neglected  [i.e.,  J  ^  =  M  +  d^u/ 

8^^)],  the  generalized  core  solutions  can  be  expressed  as 


,  8<l> 

— +  — 

and  the  current  density  components  are 


(lla) 

(llb) 


J^C  —  jqC  —  (  1  IC) 

The  methods  of  singular  perturbation  expansions  and 
matched  asymptotic  expansions  at  the  boundaries  between 
the  Hartmann  layer  and  core  region  are  applied  to  the  com¬ 
plex  system  of  coupled  equations  to  obtain  the  core  solutions 
in  the  different  core  regions  to  zeroth  order. 

For  continuity,  the  well-known  matching  principle  be¬ 
tween  the  core  and  Hartmann  boundary  layer  will  be  stated 
here  to  the  zeroth  order.  If  a  physically  realizable  function 
Uc  (x;M)  at  high  Hartmann  number  can  be  expanded  in  the 
core  as 


Uc(x;M)  =  Uco(x)  +  M  'uciix)+M  'Uc2ix)  +  ---, 

(12) 

where  x  is  an  independent  variable,  and  the  function 
(X;M)  can  be  expanded  at  iiigh  Hartmann  number  M  in 
the  Hartmann  layer  as 

u,(X-,M)  =  «6o  (X)+M-  'u, ,  {X)  -f  (X), 

(13) 

where  A'  is  a  stretched  coordinate,  then 


aco(0)  =  lim  [«,o(Ar)],  (14) 

X-  ec 

where  Uco  (0)  is  the  value  of  the  zeroth  order  function  at  the 
wall.  The  boundary  layer  value  in  the  Hartmann  layer  in 
stretched  coordinates  is  the  value  away  from  the  wall  as 
X->  oa  (see  Appendix  B  for  details). 

Cook,  Ludford,  and  Walker'*  proved  that  if  there  is  a 
jump  to  0(  1 )  in  the  velocity  u  across  the  Hartmann  layer, 
then  the  Jumps  in  j  and  <f>  are  at  most  0{M~').  The  jump 
in  ^  is  0(  A/  “  ^ )  for  an  insulator  or  a  perfect  conductor  with  a 
normal  magnetic  field,  but  it  is  0(A/  “ ' )  for  a  perfect  con¬ 
ductor  with  a  skewed  magnetic  field.  Therefore,  to  zeroth 
order,  we  shall  assume  that  the  potential  throughout  the 
boundary  is  the  same  as  the  wall  potential.  The  asymptotic 
expansions  for  the  potential  1**  ll^c  core  area  C  are  as¬ 
sumed  to  have  the  form'* 

+  (15) 

Thus,  the  matching  condition  on  the  bottom  and  top  perfect 
conductor  of  the  channel  for  the  potential  ^c2.o  in  the  core 
C2  to  zeroth  order  is 

0  =  <^c2.o  (17  =  tan  0^) 

~  —  (tan  "b  t^c2,o(^)  (16a) 

for  —  a  cos  0<|'<a  cos  0  (stationary  perfect  conductor), 
(^0  +  J")  =  l^c2,o  [V  =  sec  0  -F  (tan  6)^  ] 

~  (sec  9  -F  tan  9^)  (^)  "b  ^c2,o  (^) 

(16b) 
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{  »  acos0  -  tin6 


for  —  a  cos  0  —  sin  0<|’<a  cos  0  —  sin  0  (moving  perfect 
conductor).  The  zeroth-order  solution  for  the  potential  in 
the  Hartmann  layers  along  each  perfect  conductor  was  not 
needed  because  the  change  in  across  the  layer  is  only 
0{M  “ ' ),  and  using  the  value  of  the  potential  at  the  wall  was 
sufficient. 

Solving  these  two  simple  simultaneous  boundary  value 
equations  ( i6a)  and  (16b)  and  simple  algebraic  manipula¬ 
tion  of  the  magnetohydrodynamic  equations  (7a)-(7c)  for 


core  Cl  results  in  the  expressions 

j^ci.o  —  jifCi.Q  —  ~  cos  9{4>o  +  (17a) 

<t>c2.o  =  cos  (9((^o  -1-  ^)r]  -  sin  0(^o  +  ^)^,  ( 17b) 

“C2,0  =  (cos  9)7)  -  sin  ^(^„  -1-2^),  ( 17c) 

“c2.o  =  -  sin  9{(^a  +  i)  for  V 

=  (tan0)|'  (bottom wall),  (17d) 


“c2o  =  1  —  sin  0(^0  -f  I")  for  r)  =  sec9  -)-  tan  9^. 

(17e) 

It  is  interesting  to  note  that  to  zeroth  order  there  is  a  jump  in 
Uc2  across  each  Hartmann  layer  for  0/0. 

From  similar  arguments  the  zeroth  order  solutions  in 


core  C 1  are 

^ijci.;  '  Ac  1.0=  O'  (18a) 

^ci.o  r^ci.o  =  1>  (18b) 

and  the  zeroth-orde.  solutions  in  core  C  3  are 

Jricj.o—^’  (19a) 

^C.1.0  =  *^C3.0  ■  ■  d'  (19b) 


B.  Discussion  of  core  solutions  of  zeroth  order 

To  zeroth  order  in  the  core  C 1  the  velocity  is  ttci.o  =  1 
throughout  the  region.  Thus,  the  core  has  the  same  velocity 
throughout  the  core  region  as  the  moving  wall.  The  two 
components  of  current  density  in  core  C 1  to  zeroth  order 
j^c  ;.o  and  i.o  both  equal  zero.  Also,  the  potential  to  zeroth 
order  throughout  core  C 1  is  (pc\.o  =4‘o  + 

To  zeroth  order  in  core  Cl  the  current  density  follows 
the  external  magnetic  field  B„  and  thus  Ac  2,0  =  0  and 
Ac2.o  =  -  cos  0((^„-f  I").  The  potential  ^(77,'^)  and  u(i7,|') 
to  zeroth  order  are  functions  of  77,^  [see  Eqs.  (17b)  and 
( 1 7c )  ] .  It  IS  interesting  to  note  that  to  zeroth  order  there  is  a 
jump  in  Me  2  across  each  Hartmann  layer  at  the  top  and  bot¬ 
tom  wall  when  0/0. 

In  core  C  3  to  zeroth  order  Uc},o  equals  zero  throughout 
the  region.  The  core  C  3  flow  thus  has  no  velocity  throughout 
the  region.  Both  components  of  current  density  in  core  C3, 

jr,cxo  andAcj.o  are  zero  as  in  core  Cl.  The  potential  <Ac3.o  is 

also  zero  throughout  the  core  to  zeroth  order. 

IV.  ZEROTH-ORDER  MATHEMATICAL  SOLUTIONS  IN 
FREE-SHEAR  LAYERS 

A.  Development  of  basic  theoretical  concepts 

The  free-shear  layers'^  in  regions  /  I  and  fl  in  coordi¬ 
nates  (x,t,^)  are  shown  in  Fig.  4.  The  corner  region  of  the 
free-shear  layer  / 1  between  the  insulating  wall  and  the  per¬ 
fect  conductor  is  at  ( /  =  0,  ^  =  —  a  cos  0) .  The  free-shear 


FIG.  4.  Free-shear  layers  m  regions  f\  and  fl  in  coordinates  (tt,  t,  |)  fir 
/I,  r  =  j;  -1-  a  sin  ft  <I>  =  -f  sin  0  —  acoiO  and  for  fl,  r  =  sec  0  — 
sm  0(tan  0  —  a)  —  r) 


layer  fl  has  the  corner  region  between  the  insulating  wall 
and  the  perfect  conductor  at  (f  =  0,  ^  =  a  cos  0  —  sin  0). 

We  shall  work  out  the  fundamental  mathematical  theo¬ 
ry  for  shear  layer  /  I;  the  theory  for  shear  layer  fl  can  be 
performed  in  a  similar  manner.  To  develop  the  theory,  the 
stretched  coordinate  ^  will  be  used.  It  is  defined  as  follows; 

^  =  M  -h  a  cos  0),  ^  =  —  a  cos  0  + 

(20a) 


(20b) 


9  _  d  d 

di  di  • 

Using  Eqs.  (7a)-(7c),  the  fundamental  magnetohydrody¬ 
namic  equations  in  the  free-shear  layer  fl  for  fully  devel¬ 
oped  flow  in  stretched  coordinate  f  and  regular  coordinate  77 
are 


.  ,  I  { d^u  ..  d^u\ 


^  d^{7),^) 

077 


I  j,j.y2  0. 


(21a) 

(21b) 

(21c) 


077  0f 

The  first  terms  in  the  series  solutions  in  the  free-shear  layer 
fl  are 


Jr,iV>P  =  A/1,0  (V.f). 

(22a) 

(22  b) 

(22c) 

M  =it/'^Ay,,o(l7.^).. 

(22d) 

From  Eqs.  (22a)-(22d)  the  following  relationships  can  be 
developed  for  the  potential  function  o  (77,^) 


«/i.o(V.^’) 

yV/i.o 


ds  ’ 

~  9<l> 

077 

0Vy|(77,^) 

9r 


(2da) 

(23b) 

(23c) 
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Using  the  conservation  of  charge,  Eq.  (21c),and£qs.  {22a) 
and  (22c),  the  following  fourth-order  partial  differential 
equation  is  derived: 

dif  ~  ^  ^ 

for  —  a  sin  —  a  sin  0  -t-  sec  9  and  —  oo  <  oo  for 
free-shear  layer  /I  (see  Fig.  5). 

The  stretch  coordinate  ^  must  extend  from  —  oo  to  oo 
in  accordance  with  the  definition  of  stretched  coordinates. 
This  figure  also  shows  the  core  potentials  and  4>ci-> 
which  must  be  matched  by  the  sides  of  the  free-shear  layer. 
4> /1.0  )  must  satisfy  the  following  boundary  conditions: 

( 1 )  Boundary  condition  at  core  C2: 

CO&  d(tf>o~  a  cos  0){r}-\- a  6)  as  ^-*oo. 

(25a) 

(2)  Boundary  condition  at  core  Cl: 

'!> /1.0  -<^0  —  cos  0  as  —  00 .  (25b) 

(3)  Boundary  condition  at  moving  perfect  conductor: 

(j) y^i  0  =  (f>o  o  cos  6, 


at  1?  =  sec  0  —  a  sin  0  for  —  oo  <  f  <  oo . 

(25c) 

(4)  Boundary  condition  at  fixed  perfect: 

<^/i.o  =0  at  57  =  —  fitsin  0  for  0<f  <  00.  (25d) 

(5)  Boundary  condition  at  Hartmann  layer  between 
free-shear  layer  /  I  and  insulating  wall  (see  Fig.  6): 

{  -  a  sin  9,^)  =  sgn(sin  9)  ^  (  -  a  sin  9,^). 

Ob' 

(25e) 

Boundary  condition  ( 5 )  will  be  derived  in  the  following  sub¬ 
section. 


0(M"h 


Hattmann  layer 
H  betwaafl 
frae  shear  layer 

comer  2 
insulating 


^  comer  region 
tOtM'b  X 
which  appears 
to  free  shear 
layer  as  source 
of  electric 
current  density 


FIG.  6.  Free-shear  layer  at  region  /I. 


ers  (see  Fig.  7).  For  Hartmann  layers  adjacent  to  core  re¬ 
gions,  the  layer  matches  any  core  velocity  and  satisfies  «  =  0 
or  1,  depending  on  the  wall.  The  change  in  u  across  the  Hart¬ 
mann  layer  to  the  wall  has  a  simple  exponential  structure. 

Figures  3  and  6  show  the  Hartmann  layer  along  the  in¬ 
sulating  wall.  To  treat  the  Hartmann  layer,  we  will  return  to 
the  original  Cartesian  coordinates  (x,y,z).  Equations  ( la)- 
( Id)  for  magnetohydrodynamic  fully  developed  flow  in  the 
original  coordinates  are 

0  =  sin  9j],  -  C0&  9j\  +  M  ,  (26a) 

y  =  -  -  (sin  9)u,  (26b) 

dy 

A  =  -  -f  (cos  9)u,  (26c) 

dz 


dy  dz 


=  0, 


where 


(26d) 


B.  Boundary  condition  on  Hartmann  iayer  and  free- 
shear  layer 

For  external  magnetic  field  angles  in  the  range  of 
O<0<arctan(2fl),  the  Hartmann  layers  play  a  relatively 
passive  role  except  when  they  are  adjacent  to  free-shear  lay- 


i  =  JyiyyZ)y  + 

(26e) 

v  =  u{y,z)\, 

(26f) 

(p  =  p{y,z). 

(26g) 

The  stretching  coordinate  for  the  y  direction  in  the  Hart¬ 
mann  layer  (tangential  direction  along  insulating  wall)  is 
approximately 


I 

I 

I 


tci  *  fo  -  "cote-H 


cinS  *  aca«9 

|«-n  =  Mc9  -  atlnd 

I 

I 

»  cote  1^0  -  acotei 
I  X  (i)  -I  atine) 

I 

I 


{  ■  -acotS 
r)  X  -atine 


FIG  5.  Free-shear  layer  at  |  =  —  a  cos  ft  /I 


FIG.  7.  Hartmann  layer  H  be .  Ateen  free-shear  layer  and  insulator  at  y  =  0, 

z—  —  a. 
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Y  =  M''-y,  y  =  M-''-Y,  =  (27) 

ay  oY 

(see  Figs.  6  and  7).  The  stretching  is  needed  because  the 
Hartmann  layer  is  adjacent  to  the  free-shear  layer  where 
variables  vary  over  0{M  “  distances  through  the  layer. 
The  stretched  coordinate  normal  to  the  Hartmann  layer  is 


Z  =  M(z  +  a),  z  — 


■a  +  Af-'Z,  —  =  M—. 

dz  dZ 


The  first  terms  in  the  series  for  the  solutions  to  the  magneto¬ 
hydrodynamic  equations  are 

u  =  M''-u„(Y,Z),  <l>  =  d>HiY,Z),  (29a) 

jy=M'i'-j^„(Y,Z),  h  =  h„(Y,Z).  (29b) 

The  Hartmann  layer  must  match  the  Oi^M ''")  velocity  in 
the  free-shear  layer,  while  the  Hartmann  layer  current  den¬ 
sity  along  he  insulator  is  of  the  same  order  of  magnitude  as 
the  velocity  The  magnetohydrodynamic  equa¬ 

tions  in  the  Hartmann  layer  in  stretched  coordinates  can  be 


written  as 


0  =  sin  djy„  + 


dz-  ’ 

-  (sin^)n,„ 


— —  =  0  for  0<K<oo  and  0<Z<oo,  (30c) 

dZ 

^  +  ^  =  0.  (30d) 

dY  dZ 

At  the  insulating  wall  Z  =  0  in  the  stretched  coordinate  and 
the  velocity  u„  ( Y,Z  =  0)  —0  and  the  current  density 
(  y,Z  =  0)  =  0.  Since  the  partial  derivative  of  with 
respect  to  Z  is  zero,  <!>,/  =  ( T ).  Therefore,  the  current 

density  can  be  expressed  as 

d6„(Y) 

jy„=  -  -(sing)«„(r.Z).  (31) 

dY 

Substituting  Eq.  (31)  into  Eq.  (30a)  results  in  the  partial 
differential  equation 


d^u„{Y,Z) 


sin^  0Ufj(Y,Z)  =  (sin  6) 


d<^„{Y) 


The  solution  which  satisfies  =  0  and  which  does  not 
grow  exponentially  as  Z-* «  is 


dd>„{Y) 

u„(Y,Z)  =  — CSC0 - [1  — exp(  —  |sin0|Z)], 

dY 


where  the  term  |  sin  0 1  has  been  used  in  case  9  is  negative  in 
the  problem.  The  component  of  the  current  density  in  thej* 
direction  is  obtained  from  Eqs.  (31 )  and  (32a) 

dd>„{Y) 

jy„{Y,Z)= - !:^^exp(-|sin0|Z).  (33) 

dY 

Using  the  conservation  of  charge  condition  Eq.  (30d)  in 
stretched  coordinates 

djiH  d^(j)if{Y)  . 

- = - ^ — exp(  —  sinpZ)  (34) 

dZ  dY^^'' 

results  in  the  solution  for  ( y,Z)  after  integrating  from 
Z  =  0  to  Z  =  Z: 


Aw  =  |csc  9 


d'4>H{Y) 


[1  —  exp(  —  |sin  0  [Z)  ].(35) 


We  now  need  to  match  the  free-shear  layer  and  the  Hart¬ 
mann  layer  solutions  at  the  boundary.  The  geometric  rela¬ 
tionship  between  the  coordinates  (.x,  rj,  £)  and  (x,  y,  z)  are 
shown  on  Fig  2  and  are 

rf  =  zsin9  +ycos9,  ( 36a) 

^  =  zcos0  — _ysin  0.  (36b) 

In  the  free-shear  layer  and  adjacent  to  the  Hartmann  layer 
z=  -a  +  M-'Z,  y  =  M-''^Y,  (37a) 

Z  =  M{z  +  a),  Y^M'i'^y,  (37b) 

and 

f  =  A/'^^(|- +acos6>).  (37c) 

The  coordinate  transformations  (36a)  and  (36b)  after  sub¬ 
stituting  Eqs.  (37a)  to  (37c)  are 

17=  —  flsin^ -t-A/"''^ycos0 -l-Af"'(sin0)Z,  (38a) 

—  aco^9  —  M~''^Y%\n9  {cos  9)Z,  (38b) 

;  =  +acos9)  =  -  rsin0  +  M-'>Hcos9)Z. 

(38c) 

In  order  to  match  a  free-shear  variable,  such  as 
^/i  (V>f)>  fo  the  corresponding  Hartmann  layer  variable 
^//(y,Z),  we  substitute  the  expressions  (38a)  and  (38c) 
into  .  We  then  use  a  Taylor  series  to  obtain  the  proper 
form  for  an  asymptotic  expansion,  namely  powers  of  M 
times  coefficient  functions  of  Y  and  Z.  This  process  gives 


tpfil  —  asm9  +  M  '^^Ycos9+M  '(sin0)Z,  —  Tsin^-I-A/  ''^cos9Z] 

=  —asm9,  —  Ysm9)  +  [M  ~*'‘Ycos9  +  M~'isin9)Z] —^  { —asm9,  —  Fsin^) 

dr] 

+  M  cos  9Z  (—a  sing,  —  Tsin  0) -(•••• 


Only  the  leading  term  (  —  n  sin  0,  —  Y sin  0),  enters 
the  matching  of  the  zeroth-order  variables.  The  other  terms 
in  the  Taylor  series  expansions  of  the  zeroth  order  free-shear 
layer  variables  enter  the  matching  at  higher  orders.  Physi¬ 


cally  the  zeroth-order  free-shear  layer  solution  sees  the  Hart¬ 
mann  layer  as  a  zero  thickness  layer  at  the  horizontal  line 
1/  =  —  a  sin  d.  In  fact,  the  Hartmann  layer  has  a  small 
0{Af  ')  thickness  which  would  enter  the  higher  order 
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matching  through  the  Af  “ cos  OZd^fx  /d^  term,  and  the 
layer  lies  along  a  sloping  line,  which  would  enter  through  the 
(AT  ~  '^*Fsin  0  +  A/  " '  sin  6Z)  d<l> y,  /drj  term. 

From  the  matching  principle  between  the  free-shear  lay¬ 
er  and  the  Hartmann  layer  along  the  insulating  wall  and 
adjacent  to  the  free-shear  layer,  we  obtain 

<^//(F)  =<^y,  (  —  osin0,  —  Fsin0).  (40) 

This  expression  is  derived  in  detail  in  Appendix  C,  using  the 
velocity  instead  of  the  potential. 

From  expression  (40)  the  following  relationship  is  ob¬ 
tained: 


- ; —  =  Sin-  6 - 

dY-  dC' 


i  —asinf),  —  Tsin  0), 


X(C=-Ysmd).  (41) 

The  z  component  of  the  current  density  is  related  to  the  jg 
and components  of  the  current  density  by 

7,  =cos6»;.. -t-sin0y,,.  (42) 

To  zeroth  order  the  normal  current  density  in  the  Hartmann 
layer  j,,f  approaches  the  normal  current  density  in  the  free- 
shear  layer 

d  (  T  )  /i,  ly  -  ly  •  ^ 

7.//= — - cos6)(A/  ''-7i/i) 

aY~ 

-f  sin  6 7,,yi  (  —  a  sin  0,  —Y sin  6).  (43) 

Therefore,  using  Eq.  (41 )  the  following  relationship  holds: 
sin  9 7',,/|.o  (  —  a  sin  0,  —  K sin  9) 

=  Icsc  9  |sin-  9  (  —  a  sin  0,  —Y sin  9). 

d^- 

(44) 

The  current  coming  out  of  the  Hartmann  layer  into  the  free- 
shear  layer  equals  the  tangential  second  partial  derivative  of 
the  potential  in  the  free-shear  layer  and  is  the  fifth  boundary 
condition  discussed  in  Eq.  (25f): 

7„/i,o(  -asin^.?) 

B  ~(b 

=  sgn(sin0) - ^^  (  —  a  sin  0,  ^).  (45) 

C.  Partial  differential  equation  and  boundary  values  for 
free-shear  layer 

We  shall  assume  the  solutions  to  the  partial  differential 
equation  (24)  and  boundary  values  (25a)-(25e)  tobeinthe 
form 

(t> ) \jn  =  <t>»- a  co%  9 (46a) 
where 

/  =  7/  -t-  a  sin  0  (46b) 

and 

sin^>0  for /I,  sgn(sin0)  =  -b  1.  (46c) 

<t>(  r,^)  is  treated  as  a  perturbation  which  approaches  zero  as 
—  00,  see  boundary  condition  (25b).  Also  <l>(t,f)  ap¬ 
proaches  zero  at  the  moving  perfect  conductor  [see  condi¬ 
tion  (25c)  !•  The  partial  differential  in  terms  of  •Pir.p  with 
its  five  boundary  conditions  in  terms  of  /,  f  is 


. . . 

— -  ,  =  — — T —  for  0</<sec  9,—  co  <g<oo; 

at-  ay 

(47a) 

Boundary  condition  at  core  C2: 

*P(t,g)  -*  (^0  -  a  cos  9)  (t  cos  9  —  1),  as  f-oo; 

(47b) 

Boundary  condition  at  core  C 1 : 

<l>(f,f)—0,  as  g—  —  co;  (47c) 

Boundary  condition  at  moving  perfect  conductor: 

<^(t,g)=0,  at  r  =  sec0,  for  —  oo  <g<  oo;  (47d) 
Boundary  condition  at  fixed  perfect  conductor: 

=  —  ((^n  —  a  cos  9),  at  7  =  0  for  0<^<  oo; 

(47e) 

Boundary  condition  at  Hartmann  layer  adjacent  to  free- 
shear  layer: 

(47f) 

We  solved  the  partial  differential  equation  and  its 
boundary  values  in  terms  of  the  modified  complex  Fourier 
transform in  the  form 


<J>(a) 


exp(  -  /Ag)<i>(t,g)dg, 


where  k  is  generalized  to  become  complex  in  the  lower  half¬ 
plane  (i.e.,  k=:/c,  —  ik^),  or  Im(/:)  <0.  Accordingly,  the 
modified  Fourier  inversion  formula"*  is  written  as 

<l)(r,f)=-^r  ‘  t%^{ikgmtM)dk.  (49) 

(277)  7-«-n 

The  inversion  contour  is  in  the  lower  half-plane,  along  the 
line  Im(^)  =  —  r. 

The  Fourier  transform  of  the  partial  differential  equa¬ 
tion  is  in  the  form 

—r  —  k*^,  for  O<7<sec0,  (50) 

Bt " 

where  integration  by  parts  has  been  used  to  obtain  the  above 
form.  The  right-hand  side  requires  four  integrations  by  parts 
with 


U\^[-ikg)^(r],g)\  =0. 

\  wC  /  f 


The  upper  limit  f  —  oo  gives  zero  because  Im(^)  <0;  the 
lower  limit  f  —  —  oc  gives  zero  because  of  condition  (47c). 
_  Fourier  transforming  Eq.  (47d),  we  recognize  that 
<l>  =  0,  at  7  =  sec  0,  and  assume  that  the  general  solution  of 
the  Fourier  transformed  partial  differential  equation  [e.g., 
Eq.  (50) )  is 

‘l>(7,^)  =/l(A')sinh|^-(sec0-  7)],  (52a) 


where 


^(7,^) 


=  —  A:',4(A:)cosh[A:-(sec  0  — 7)1  (52b) 


d-^{t,k) 


=  A‘‘/4(A)sinh|/:  ’(sec  9  —  t)]. 
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The  boundary  value  of  Eqs.  (47f)  can  be  written  in  the  more 
general  form 


d^-  at 


where 


°°  <^<0> 

l/(f),  for  0<f<oo, 


{53a) 


(53b) 


and  /  (f )  is  an  unknown  function.  The  Fourier  transform  of 
F  (f )  is  expressed  as 


F{k)=~[,0,k)-k^-<^m). 

dt 


(54) 


Substituting  the  general  expression  for  d>(/,/:),  Eq.  (52a), 
results  in  the  expression 

F  {k)  ~  —  k‘^(k)[cosh(/:- sec  0)  +sinh(k'sec0)], 

(55a) 

where 

A{k)  =  —  F  (/c)/k '  exp{/c  ’  sec  0).  (55b) 

Therefore  ^(t,k)  can  be  conveniently  expressed  as 

<I>(r,k)  =  —  F(/L)sinh[k“(sec  0  — /)]/k^exp(k-sec  0) 

(56a) 

=  -  [  F(k)/2k"]{exp(  -rk') 

—  exp(  —  k  ^(2  sec  <9  — /)  ]},  (56b) 

where  F(/:)is  an  unknown  function. 

The  convolution  theorem  can  be  used  to  obtain  d>(r,f ) 


<!>( 


/.f)  =  r  "  F(^*)(?(t,f-f‘)£(r  (57a) 

J  v.  ■  tf 


and 


(57b) 


using  the  definition  of  F  (p  expressed  in  Eq.  (53b). 

The  Green’s  function  G(t,p  can  be  obtained  by  the 
modified  inverse  Fourier  transform  [see  Eqs.  (57a)  and 
(57b)]; 

G(t,p=---j—\  {exp(-rk-) 

(4ir)  J  -  *  r, 

dk 

-  expl  -  k-{2  sec  0  —  f)  ]}exp(/kf)  — - . 

k  - 
(58) 


There  are  no  poles  along  the  real  k  axis  at  k  =  0,  which  can 
be  proven  by  expanding  the  exponentials  in  power  series  and 
noting  that  these  terms  cancel  the  k '  term  in  the  denomina¬ 
tor.  Therefore,  the  integration  of  Eq.  (58)  can  be  performed 
with  real  k.  Expanding  the  exponential  (i.e., 

=  cos  k^  +  i  sin  kp  the  integration  can  be  changed  from 
—  00  to  00  to  0  to  00  by  neglecting  the  odd  function  or  sine 
term. 

G(tp  was  integrated  by  parts  to  obtain  an  integral 
form  which  can  be  simply  integrated  to  exponential  and  er¬ 
ror  functions 


G(t,p=—\  {/exp(  —  rA:^)  —  (2sec0— 0 

IT  Jo 

Xexp[  —  k^(2  sec  $  —  t)]}cos  dk 

+  — ^  f  {(exp(-rA:^) 

(.2ir)  Jo 

—  expl  —  A: ^(2  sec  0  —  /) ]}  dk. 

(58') 

It  should  ')e  noted  that  sin  k^  /k  does  not  produce  a  singular¬ 
ity  ?Lsk-*0,  and  that  integration  along  real  k  is  still  allowa¬ 
ble.  G(t,p,  Eq.  (58),  was  evaluated  by  using  the  two  inte¬ 
grals 

Jr»oc 

exp(  —  tk^)cos(kpdk 

0 

=  (7r''-/2)/-''*exp( -fV4/),  (58a) 

^ec 

exp(  —  /A:  ^)  (sin  k^ /k)dk 
Jo 

=  (7r/2)erf[if/-'^=]  (58b) 

to  have  the  form 

G(l,p 

=  ^v-'^Ht'^^expl  ->pt~']  -  i2sec9-t)''‘ 

X  expl  -  if^(2sec0-/)-']} 

-f  (f/4){erf[^f/-'«] 

-erf[Jf(2sec0-/)-'^2j;} 


Remembering  that  d>(/,f )  can  be  expressed  by  the  con¬ 
volution  integral  [see  Eq.  (57b)] 

f(^*)G(U-^*)d^\  (59) 

Jo 

We  must  also  satisfy  the  boundary  condition  at  the  fixed 
conductor 

•I*  =  -  (^0  —  0  cos  0),  at  f  =  0  for  0<f  <  00.  (60) 
Therefore,  the  following  relationship”  must  hold. 

Jo 

=  -  (<^()-acos0)  for  0<f<oo.  (61) 

Now  all  the  boundary  conditions  have  been  accounted  for  in 
the  problem. 

In  Eq.  (61 ),  we  will  change  coordinates  to  reduce  the 
integral  equation  to  a  form  that  is  independent  of  sec  dmd  of 
(^„  —  a  cos  9)..  The  following  transformation  will  be  per¬ 
formed: 


^  =  2{2sec9)''-Z,  =  2{2sec9)''^Z*,  (62a) 

d^*  =  2{2  sec  9)''UZ*,  where  0<^-,  Z<<»,  (62b) 
/[2(2sec0)''’-Z] 

=  fi0=  -[(<t>o-acos9)/2sec9]F(Z).  (62c) 
The  integral  equations  is  now  in  the  form 
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F{Z*){\Z-Z*\tdc{\Z~Z*\) 

Jo 

-7r'/‘exp[-(Z-Z*)-]VZ*  =  l  for  0<Z<oo. 

(63) 

This  equation  must  be  solved  for  f  ( Z)  in  order  to  solve  for 
the  velocity  profile  directly  in  the  free-shear  layer  (see,  for 
example,  Eqs.  (22d)  and  (23a)  (.This  work  will  be  reported 
in  the  third  part  to  this  paper,  which  is  in  progress. 


D.  Characterization  of  free-shear  layers  f  1  and  /2 

For  the  free-shear  layer  /I  for  O<0<arctan(2a),  we 
will  use  a  coordinate  system  (x,  r,  J")  (see  Fig.  4)  with  the 
origin  at  >>  =  0  and  z  =  —  a,  with  ^  perpendicular  to  B  and  t 
parallel  to  B.  The  free-shear  layer  / 1  intersects  the  surface 
of  the  moving  perfect  conductor  at  =  1  and  2=  —a 
-{-  tan  0,  where  the  conductor’s  electric  potential  is 
*1^/1  =  <^()  —  ^cos  0  =  sin  0  — acos  0.  Here  is 
the  moving  conductor’s  electric  potential  at  the  middle  of  its 
surface,  i.e.,  at  =  1  and  z  =  0.  The  potential  of  the  mov¬ 
ing  conductor’s  surface  is  then  (l>  =  <pm  +  z  cos  6  =  <f>„ 
-f  sin  0  -i-  at  >>  =  1.  If  =  0,  the  potential  difference 
between  the  moving  and  fixed  electrodes  is  zero  at  z  =  0, 
while  the  potential  of  the  moving  electrode  is  greater  than  or 
less  than  that  of  the  fixed  conductor  for  z  >  0  or  z  <  0,  respec¬ 
tively.  On  the  left-hand  side  of  free-shear  layer  /  I,  the  po¬ 
tential  in  core  Cl  is  for  all  t.  On  the  right-hand  side  of 
/I,  the  potential  in  core  C2  varies  linearly  from  <l>y,  at 
t  =  sec  6  to  zero  at  r  =  0.  Therefore,  there  is  a  jump  in  poten¬ 
tial  across  / 1  and  the  magnitude  of  thisjump  varies  linearly 
from  zero  at  r  =  sec  0  to  —  <l>y ,  at  r  =  0.  Since  the  thickness 
of /I  isO(3/“''’),(?^/<?|must  be  large,  namely  0(A/''‘), 
inside /I.  In  the  |  component  of  ohms  law. 


d<p 


+  u. 


(64) 


4  is 0(M  [see Eq.  (22c)  ],  so u  =  d^/d^  =  0{M ''*). 

This  leads  to  the  conclusion  that  the  velocity  in  the  free- 
shear  layers  is  much  larger  than  the  velocity  of  the  moving 
wall.  Therefore,  even  though  the  boundary  layer  is  thin 
0(M  it  involves  large  velocities,  so  that  the  total  flow 
inside  the  free-shear  layers  is  comparable  to  that  in  cores  C 1 
and  C2: 


I, 


«  dj's  —  (1  — /cos  0).  (65) 

/across/I 

From  this  equation  it  can  be  deduced  that  if  ,  >0,  the 
free-shear  layer  / 1  is  a  high  velocity  jet  in  the  minus  x  direc¬ 
tion.  Conversely,  if  <!>/,  <0,  /I  is  a  jet  in  the  plus  x  direc¬ 
tion.  Noting  Eq.  (65),  it  is  interesting  to  note  that  the  inte¬ 
gral  varies  linearly  from  0  at  the  moving  electrode  to  a 
maximum  of  — ‘b/i  at  the  comer  at  z  =  —a,y  =  0.Qf,, 
the  total  flow  inside /I  inplusxdirection  (i.e.,<t>y-,  <0),is 


0  r 

Qf^  =  \  u  dt^  -  sec  e.  (66) 

Jo  Jacross/i 

Also,  for  the  discussion  of  free-shear  layer  /2,  we  will 
use  a  similar  coordinate  system  t,  x)  (see  Fig.  4).  We 
shall  define  ^  f  2.  =<f>m  +  a  cos  0  as  the  potential  at  the  in¬ 
tersection  of  /2  and  the  surface  of  the  moving  electrode,  i.e.. 


My  =  1  and  z  =  a.  By  the  previous  arguments  for  / 1  it  can 
be  shown  that  the  velocity  in  the  free-shear  layer  is  much 
larger  than  the  velocity  of  the  moving  wall.  The  integral  of 
the  velocity  across  the  boundary  layer  is 


I 


'ACROSS/2 


ud^=  —  t<t>f2  COS  0. 


(67) 


If  ^f2  >  /2  is  a  high  velocity  jet  in  the  minus  x  direction. 

If  <l>y2  <  0,  /2  is  a  jet  in  the  plus  x  direction. 

The  integral  [Eq.  (67)]  varies  linearly  from  0  at  the 
fixed  electrode  to  a  maximum  of  —  <l>y^2  Ihc  corner;'  =  1, 
2  =  a.  Qf2<  the  total  flow  in  the  plus  x  direction  (i.e., 
<l>^2  <0),is 


in 
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u  rf^d/s  —  (<l>o  sec  0.  (68) 


'  ACROSS /2 


The  dimensional  flow  is  determined  by  the  definition 


Q*  ^  (.U,dL'-)Q. 


(69) 


The  total  dimensionless  flow  in  the  x  direction  equals 
the  total  flow  in  the  three  core  areas  Qc  1 ,  Qci  >  and  Qa  and 
the  total  flow  in  the  two  free-shear  layers  Qfi  and  Qf2  ■  The 
flows  in  the  Hartmann  layers  of  0(Af  “' )  will  be  neglected. 
Therefore,  the  total  flow  Q  is  expressed  as 

Q=Qct  +  Qc2  +  Qci  +2/1  +  2/2  +  0{M~') 


-f  Hartmann  layers,  (70a) 

where 

2/1  =  -('J>/.(sec0),  (70b) 

2/2  =  -  (<l>/2  (sec  0),  (70c) 

eci=Uan6),  (70d) 

Qc2  =  (2acos0  — sin0)(4cos0  — tan^),  (70e) 
2c3=0.  (700 

After  substitution  of  these  expressions  and  simplification, 
the  total  flow  in  the  x  direction  Q  is 

Q  =  cos  0(.a  cos  0  —  ( sin  6)  —(f)„{2a  sin  9  -i-  cos  6) 

-f  0(4/“')  +  Hartmann  layer  terms.  (71) 


Let  us  now  examine  the  properties  of  the  velocity  pro¬ 
files  for  no  net  current  flow  between  the  electrodes.  The 
0(1)  electric  currents  in  cores  C 1  and  C  3  are  blocked  by  the 
insulating  sides,  so  that  the  “active”  part  of  the  moving  elec¬ 
trode’s  surface  is  from  2=  —  a  -1-  tan  0  to  z  =  a.  For  zero 
net  current,  the  potential  at  the  middle  of  this  active  part  of 
the  moving  electrode  must  be  zero,  so  that 

<i„=-(sin0.  (72) 

Then  =  —  (a  cos  0  j  sin  0)  which  is  negative  be¬ 
cause  tan  9  <  2a.  Therefore,  / 1  is  a  jet  in  the  plus  x  direc¬ 
tion.  Also,  ^f2  =  (acos0  — Jsin0)  =  —  >0,  and 

thus  /2  IS  an  equal  jet  in  the  minus  x  direction.  The  total  Q  is 
equal  to  a  [seeEq.  (70a)]. 

At  this  point,  it  is  interesting  to  note  for  any  </>„ ,  Uc2 
varies  linearly  with  z  from 


«C2  =y  —  ^n  (sin  ^)>  at  z=  —  a  -f^'tan  9  (73) 

to 

Uc2  —  y  —  '^/2  (sin  0),  at  z  =  a  —  (1  —  >')tan  9.  (74) 
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Now  we  shall  plot  a  typical  plot  of  Uc2  vs  z  for  a  fixed 
value  of  y  for  zero  net  current  ( Fig.  8 ) .  The  jets  in  the  free- 
shear  layers  in  / 1  and  /2  exactly  cancel  each  other  [  see  Eq. 
(66)].  Note  that  Uc2  at  position  1  can  be  greater  or  less  than 
1, depending on>'and0(i.e.,Mc2  >;').  At  position  2,  ^^2  can 

be  greater  or  less  than  0,  depending  on;'  and  0(i.e.,«c2>>')- 

Let  us  now  evaluate  the  characteristics  of  the  velocity 
profile  in  terms  of  quantity  of  flow  for  a  net  transport  current 
from  the  moving  perfectly  conducting  electrode  to  the  per¬ 
fectly  conducting  stationary  electrode.  The  net  transport 
current  is 

/=  (Ax)  J  jy  dz<0,  for  <{>„>  -k  sin  6.  (75) 

<i>„  is  the  mean  voltage  at  the  moving  electrode.  As  in¬ 
creases  from  —  1  sin  0  the  jet  in  the  positive  x  direction  in¬ 
side  the  region  /I  (i.e.,  +  sin  ^  —  acos<0)  de¬ 

creases  in  magnitude  (the  quantity  of  flow  Qj  ^  decreases  in 
magnitude).  All  velocities  in  core  region  C 2  decrease.  Also, 
the  jet  in  the  negative  x  direction  inside  region  /2  (i.e., 
^  cos  0>O)  increases  in  magnitude.  When<^„ 
reaches  the  value  (a  cos  0  —  sin  6),  then  ^ =  0,  so  that 
there  is  no  jet  inside  region  /l,and«c2  =J',  at2  =  —a+y 
tan  6.  At  this  .  ’i^/2  =  2a  cos  ^  —  sin  6;  so  that  the  flow 
in  the  negative  x  direction  inside  region  /2  is  two  times  what 
it  was  for  =  —  (  sin  6.  The  total  flow  Q  equals  a  —  a  cos 
(20)  —  (a’ —  i)sin  2ft  Finally,  as  d)^  increases  further 
from  a  cos  0  —  sin  0,  jets  in  both  regions  / 1  and  fl  are  in 
the  negative  x  direction. 

Let  us  now  investigate  the  velocity  profiles  for  a  net 
transport  current  from  the  stationary  electrode  to  the  mov¬ 
ing  electrode  as  expressed  as 

/=(Ax)J  y^£(2>0,  for  < -Jsinft  (76) 

First,  d>m  decreases  from  -  ( sin  0  to  larger  negative 
values,  with  the  result  that  the  jet  in  region  / 1  increases  in 
volume  of  total  flow.  All  velocities  in  region  C  2  increase. 
Also,  the  jet  in  the  negative  x  direction  in  region  /2  de¬ 
creases  in  magnitude.  When  <f>„  =  —  a  cos  0,  so  that 
(py,  =  —  (2acos  0  — sin  0),sothattheflowinthepositive 
X  direction  inside  /  I  is  two  times  that  for  (f)„  =  sin  0. 
For  this  ij>„,  d>y2  =0,  so  there  is  no  jet  inside  fl  and 
Uci  =  >>,  at  z  =  a  —  ( 1  —  y)  tan  0.  The  total  flow 

Q  =  a-\-  a  cos(20)  -t-  (a^  —  j)sin  16. 


FIG  8  Plot  of  typical  1/ vs  2  for  fixed  value  of  >>  for  zero  net  current. 
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Finally,  as  d>m  decreases  further  from  —  a  cos  0,  both  f\ 
and  fl  jets  are  in  the  positive  x  direction. 

From  Eqs.  (70a)  and  (70b),  the  difference  in  total  vol¬ 
ume  of  flow  in  regions  f\  and  fl  is 

Qf^  -  Q/2  ■■■=  U2a  -  tan  0)  >0  (77) 

so  that  the  difference  is  independent  of  6„ . 

V.  EXPERIMENTAL  EVIDENCE  FOR  FREE-SHEAR 
LAYERS 

Since  we  have  analyzed  the  free-shear  layers  mathemat¬ 
ically  by  fundamental  magnetohydrodynamic  theory  and 
characterized  some  of  the  basic  properties  of  these  layers,  we 
must  ask  whether  there  is  any  experimental  evidence  of  jets 
of  boundary  layer  thickness  0{M  "  The  answer  is  yes. 
Hunt  and  Malcolm”  ''*  treat  theoretically  in  cylindrical  co¬ 
ordinates  and  experimentally  the  problem  of  a  circular  elec¬ 
trode  embedded  in  an  insulating  wall  in  „  high  external  mag¬ 
netic  field  in  a  stationary  fluid.  They  predict  a  free-shear 
layer  as  a  jet  in  the  azimuthal  direction  between  the  core 
region  and  the  electrode,  and  between  the  core  region  and 
the  insulating  wall.  They  predict  that  Vg  =  The 

free-shear  layer  is  a  jet  in  the  0  direction,  but  the  jet  is  in  the 
-f  0  direction  for  z>0  and  the  —  0  direction  for  z<0. 
(z  =  0  is  midway  between  the  two  circular  electrodes. )  In 
their  work  experimental  measurements''*  show  evidence  of 
these  jets. 

However,  that  work  is  somewhat  different  from  the  jets 
we  discuss  here,  in  which  each  jet  must  be  entirely  in  the  -t-  x 
or  —  X  direction.  For  the  work  herein,  the  jet  goes  from  zero 
at  the  intersection  with  the  perfect  conductor  to  a  maximum 
at  the  corner  where  the  insulator  and  other  conductor  meet. 
Both  cases  involve  linear  variation  of  velocity  along  the 
length  of  the  layer. 

VI.  TOTAL  CURRENT  AND  POWER  LOSSES  BETWEEN 
ELECTRODES  IN  THE  CORES 

A.  Development  of  basic  theory 

We  derived  the  total  current  between  the  electrodes  (see 
Figs.  1  and  2)  for  the  core  regions.  In  the  core  regions  Cl  and 
C  3  the  current  densities j^c  i  i  Jijcj  >  j^ci  were  deter¬ 
mined  to  be  zero  to  0(  1 ) ,  when  neglecting  terms  oiO{M~') 
in  the  magnetohydrodynamic  equations.  In  the  free-shear 
layers  /land  /2,  the  current  densities;',/,  andy,,y-2  were 
determined  to  be  of  (7(1).  But  since  the  thickness  of  the 
boundary  layer  is  0{M  “  '^‘)  the  total  current  is  0(M  "  ''^) 
in  the  i/  direction.  The  current  densities  in  the  core  C  2  were 
derived  to  be  j^ci  =  0>  and  y,c2  =  —  cos  0(<^o  + 1”)  to 
(7(1),  neglecting  terms  lo  0{M  Here,  (f>g  is  the  nondi- 
mensional  potential  of  the  perfectly  conducting  moving  wall 
at  ^  =  0,  in  reference  to  the  stationary  perfectly  conducting 
wall  at  the  bottom  of  the  channel. 

The  zeroth  order  total  current  /  to  (7(  1 )  carried  by  the 
core  C2  is 

Jacos  6—s,\n0 

j„c2  (78) 

-  ocos  0 

where  Ax  =  nondimensional  axial  length  along  channel, 
/  =  nondimensional  total  current  in  Cl,  1',=  —  a  cos  6 
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(coordinate  for  comer  2),  and  I",  =  a  cos  0  —  sin  0  (coordi¬ 
nate  for  comer  3).  The  total  current  terms  neglected  in  the 
two  free-shear  layers  and  the  two  Hartmann  layers  along  the 
insulatingwallsareofO(3/  “  '‘-)m60(M  respectively. 

Integration  of  the  expression  for  the  total  nondimen- 
sional  current  I  in  core  C2  is 

7  =  —  ( Ax)cos  +  4  sin  6)  (la  cos  6  —  sin  6), 

(79) 

where  is  the  mean  voltage  of  the  moving  wall  and  is 
defined  as  —  sin  =d>3ity  -1  and  z  =  0,  the 

potential  at  the  center  of  the  moving  wall.  The  new  dimen¬ 
sional  quantities  with  superscripts  (♦)  are  defined  as 

=  (80) 
=  (81) 

The  other  variables  were  defined  previously  in  the  paper. 

The  total  nondimensional  zeroth-order  total  Joulean 
power  loss  Pj  between  the  electrodes  to  0(Af ')  in  core  Cl 
can  be  derived  by  integrating  the  current  density  squared  in 
the  7]  direction,  in  core  region  Cl: 

/•a  COS  ~  d 

P, -(Ax)M-(sec0)  (82) 

J  a  cos  0 

where  the  sec  0  comes  from  the  integral  with  respect  to  rj. 
The  total  zeroth  order  Joulean  power  loss  terms  neglected  in 
the  two  free-shear  layers  and  two  Hartmann  layers  along  the 
insulating  walls  are  0(M^'-)  and  0(M),  respectively. 

The  expression  for  the  power  loss  Pj-  after  integration  is 

P,  =  (Ax)A/-  cos  d(la  cos  6  —  sin  6)  [<4,„(<i„  -I-  sin  6) 

+  J(a' cos’ P  — a  cos  0sin  0  -f  sin’P)],  (83a) 
where  the  dimensional  Joulean  power  loss  is  defined  as 
P1  =  {(i;Ul)P,.  (83b) 

The  current  densities  in  the  cores  C 1  and  C  3  were  deter¬ 
mined  to  be  zero  to  0(  1 )  and  thus  do  not  contribute  to  the 
zeroth  order  total  current  and  zeroth  order  Joulean  power 
losses. 

It  should  be  noted  that  the  expression  for  P  f  when  the 
magnetic  field  becomes  transverse  P  =  0,  and  no  external 
potential  is  applied  (fi„  =  0;  approaches  the  result  obtained 
by  the  authors  when  working  the  corresponding  problem 
using  the  infinite  series  type  of  mathematical  solutions  pre¬ 
sented  by  Hughes  and  Young'" 

P,  =\(^x)-M'a\  (84) 

The  nondimensional  resistance  of  the  rectangular  chan¬ 
nel  per  unit  length  R  in  core  region  C  2  is  defined  as  the 
nondimensional  length  of  the  resistor  L  divided  by  the  non- 
dimensional  cross-sectional  area  A  of  the  resistor 

R^L/A  =  \/(la  -  tan  6')(cos  ■e)(^x),  (85) 

where  A  =  (la  —  tan  6)cos  6(Ax)  and  L  =  sec  6.  The  di¬ 
mensional  resistance  /?  *  is  defined  as 

R*^(\/La)R.  (86) 

B.  Results  and  discussion 

Figures  9-14  present  straight  line  graphical  data  for  the 
nondimensional  total  current  per  unit  channel  length  in  the 
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FIG  9.  Total  current  per  unit  length  of  ciiannel  vs  positive  mean  potential 
for  aspect  ratio  a  =  \. 

core  Cl,  (I /Ax ) ,  versus  the  nondimensional  mean  voltage  at 
the  center  of  the  moving  conductor  at  the  top  of  the  rectan¬ 
gular  channel  (f>„ .  The  potential  is  referenced  to  the  per¬ 
fect  conducting  wall  at  the  bottom  of  the  rectangular  chan¬ 
nel,  which  is  assumed  to  have  the  potential  <^  =  Q.  Equation 

(79)  is  used  to  numerically  calculate  the  data.  Equations 

(80)  and  (81 )  show  the  parameters  used  to  dimensionalize 
the  nondimensionalvariables. 

Figures  9-1 1  present  straight  line  plots  of  positive  mean 
voltage  <{>„  plotted  versus  total  current  per  unit  length  (// 


FIG.  10.  Total  current  per  unit  length  of  channel  vs  positive  mean  potential 
(#,„  for  aspect  ratio  a  =  2. 
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FIG.  1 1  Tolal  current  per  unit  length  of  channel  vs  positive  mean  potential 
<},„  for  aspect  ratio  a  =  5 

Ax).  For  the  range  of  variables  studied  herein  the  (I /Ax) 
values  are  all  negative  quantities.  The  straight  line  plots 
show  that  for  positive  mean  voltages  the  total  current  is  flow¬ 
ing  from  the  moving  wall  to  the  bottom  perfect  conductor  in 
the  channel  in  core  C  2.  Each  figure  presents  data  for  a  partic¬ 
ular  aspect  ratio  ( i.e.,  a=  1 , 2,  5 ) .  In  each  figure  the  angle  6 
IS  placed  along  each  straight  line  plot  where  6  shows  the 
external  magnetic  field  angle.  The  straight  line  plots  show 
that  as  the  positive  mean  voltage  <!>,„  increases  and  the  exter¬ 
nal  magnetic  field  angle  6  decreases,  the  total  negative  cur¬ 
rent  per  unit  length  in  the  rectangular  channel  increases  for 


=  mean  potential 


FIG,  12  Total  current  per  unit  length  of  channel  vs  negative  mean  potential 
d,,,  for  aspect  latio  a  =  1 
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FIG  13.  Total  current  per  unit  length  ofchannel  vs  negative  mean  potential 
d„,  for  aspect  ratio  a  =  2. 

the  range  of  variables  interest.  The  total  current  flows  from 
the  top  moving  electrode  to  the  bottom  fixed  electrode. 

Figures  12-14  presents  straight  line  plots  of  negative 
voltagesi^,,,  plotted  against  total  current  pair  unit  length  (// 
Ax)  for  core  C  2.  For  the  range  of  variable  of  interest  herein 
the  (//Ax)  values  are  all  positive  quantities.  These  straight 
line  plots  show  that  for  negative  mean  voltages  the  total  cur¬ 
rent  is  moving  from  the  bottom  stationary  perfect  conduct¬ 
ing  wall  to  the  moving  perfect  conductor  at  the  top  of  the 
channel  in  core  C  2.  Each  figure  presents  data  for  a  particular 
channel  with  a  given  aspect  ratio  (i.e.,  a  =  1,  2,  5).  In  each 
figure  the  external  field  angle  9  is  placed  alongside  each 


^  mean  potantial 


FIG  14.  Total  current  per  unit  length  of  channel  vs  negative  mean  potential 
for  aspect  ratio  a  =  5 
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straight  line.  These  plots  show  that  as  the  negative  mean 
voltage  <f>„  increases,  as  the  aspect  ratio  of  the  rectangular 
duct  a  increases,  and  as  the  external  magnetic  field  angle  6 
decreases,  the  total  positive  total  current  per  unit  length  of 
the  duct  increases  for  the  range  of  variables  of  interest. 

It  is  interesting  to  note  at  this  point  that  positive  mean 
potentials  <!>„  give  a  larger  absolute  value  total  current  per 
unit  length  values  j  { I  /  Ax )  |  than  the  corresponding  negative 
potential  in  a  similar  channel  with  the  same  6>>0.  Each 
straight  line  intersects  the  /  =  0  axis  =  —  i  sin  6. 

In  Fig.  15,  curves  of  nondimensional  resistance  of  the 
duct  per  unit  length  R  in  core  C  2  are  plotted  versus  external 
magnetic  field  angle  9.  Three  curves  are  shown  with  aspect 
ratios  of  a  =  1,  2,  and  5.  Each  curve  has  the  lowest  total 
resistance  at  the  smallest  value  of  9  and  increases  as  9  in¬ 
creases.  The  resistance  of  the  duct  decreases  with  aspect  ra¬ 
tio  a,  as  would  be  expected.  The  nondimensional  resistance 
of  the  rectangular  channel  per  unit  length  in  core  C2  is  de¬ 
fined  as  the  resistor  length  L  divided  by  the  cross-sectional 
area  A  of  the  resistor  and  is  derived  in  Eq.  ( 85 ) .  The  dimen¬ 
sional  resistance  R  per  unit  length  is  defined  in  Eq.  (86). 

Figures  16  to  19  show  the  nondimensional  Joulean  pow¬ 
er  losses  per  unit  length  of  channel  Pr  versus  the  mean  vol¬ 
tage  at  the  center  of  the  moving  perfectly  conducting 
wall  (i.e.,  moving  electrode)  for  different  external  magnetic 
field  angles  9.  For  both  positive  and  negative  potentials  the 
power  losses  are  positive,  since  heat  is  lost  from  the  system 
due  to  Joulean  heating.  Power  losses  for  external  magnetic 
field  where  9  equals  zero  (i.e.,  transverse  external  magnetic 
field)  were  reported  in  an  earlier  paper  by  Brown  et  al."  It 
should  be  noted  that  the  power  loss  value  from  the  negative 
potential  (p„  is  slightly  less  than  the  analogous  power  loss 
value  from  the  corresponding  positive  potential.  This  is  as 
expected  because  the  positive  potential  drives  a  larger  abso¬ 
lute  value  of  the  total  current  per  unit  length  of  the  channel 


$  IdcgrMt) 


FIG.  15.  Nondimensional  resistance  per  unit  length  of  channel  vs  angle  of 
external,  homogeneous  magnetic  field  6. 
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FIG.  16  Nondimensional  Joulean  power  losses  per  unit  length  of  channel 
versus  positive  mean  potential  ^ „  for  aspect  ratios  o  =  1.2,  and  5  and  exter¬ 
nal  magnetic  field  0  =  10",  30'.  and  60* 

than  the  analogous  negative  potential.  (See  plots  of  the  total 
currents  per  unit  length  of  the  channel  in  the  preceding  part 
of  this  section. ) 

VII.  DISCUSSION  AND  CONCLUSIONS 

Fully  developed  viscous  liquid-metal  flows  in  a  rectan¬ 
gular  channel  (duct  flow)  with  insulating  side  walls  that  are 
skewed  to  a  high,  homogeneous  external  magnetic  field  were 
studied.  A  perfectly  conducting  moving  top  wall  with  an 


=  mean  potantlal 


FIG  "17.  Nondimensional  Joulean  power  losses  per  unit  length  of  channel 
vs  negative  potential  for  aspect  ratios  a  =  1,2,  and  5  and  external  mag¬ 
netic  field  angles  6  =  10",  30",  and  60*. 
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FIG.  1 8.  Wall  w  ilh  the  unit  normal  ij  at  an  angle  0  to  the  e.xternal  magnetic 
field  direction 


external  potential  and  a  stationary  perfectly  conducting  bot¬ 
tom  wall  at  zero  potential  act  as  electrodes  and  are  also 
skewed  to  the  external  magnetic  field.  The  angles  of  the  ex¬ 
ternal  magnetic  field  B,,  is  in  the  range  of  values  0<6<  arc- 
tan  la.  The  flow  in  the  channel  is  characterized  by  high 
Hartmann  number  M,  high  interaction  parameter  N,  low 
magnetic  Reynolds  number  R,„ ,  channel  aspect  ratio  a,  and 
external  magnetic  field  angle  d. 

In  this  work,  since  the  Hartmann  number  is  always  as¬ 
sumed  large  (i.e.,  1 ),  it  was  determined  by  magnetohy¬ 

drodynamic  theory  that  the  channel  could  be  divided  into 
subregions  shown  in  Fig.  3.  In  the  three  channel  core  subre¬ 
gions  C  1 ,  C  2,  and  C  3  the  first  partial  derivatives!?  /<?|;and  <?  / 
drj  are  of  0{  1 ).  The  six  Hartmann  layers,  h  \-h  6,  along  the 
sides  of  the  duct  are  of  OiM  ~ ' )  thickness  and  have  0(A/) 
normal  derivatives.  /  I  and  fl  are  two  free-shear  layers  (or 
interior  regions)  which  have  unique  characteristics.  The 
free-shear  layers  separate  the  three  core  regions  from  each 
other.  A  free-shear  layer  lies  along  the  magnetic  field  lines 
through  the  appropriate  channel  corner  and  is  driven  by  a 
current  density  singularity  at  this  corner  (see  Figs.  4-6). 
The  layers  have  thickness.  Theflrst  partial  de¬ 

rivatives  are  of  order  <?/(?^  =  0{M''-)  m6d/d'q  =  0{\) 

The  general  mathematical  solutions  in  the  cores  were 
obtained  from  the  basic  magnetohydrodynamic  equations 
( 7b )  and  ( 7c )  in  the  core  if  the  term  of  0(  A/  “  ’ )  is  neglected. 
The  method  of  singular  perturbation  expansion  and 
matched  asymptotic  expansions  at  the  boundaries  between 
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the  Hartmann  layer  and  core  region  is  applied  to  the  com¬ 
plex  system  of  coupled  equations  to  obtain  the  core  solutions 
in  the  different  core  regions  to  zeroth  order.  It  is  assumed 
that  the  physically  realizable  functions  in  the  core  Uc  (x\M) 
can  be  expanded  at  high  Hartmann  number  3/ as 

Ucix;Af)  =  Uco(x)  +  (l/A/)Uc-i  (x) 

+  (1/M-)Uc2(x)+-,  (87) 

where  x  is  an  independent  variable.  The  function  {x;M) 
can  be  expanded  at  high  Hartmann  number  M  in  the  Hart¬ 
mann  layer  as 

(X\M)  =  Kfco  (■*)  +  ( 1/Af  )«i  1  (.X) 

+  (1/M-)u,2(X)+  -,  (88) 

where  X  is  a  stretched  coordinate.  The  zeroth  order  terms 
are  I/co  ^nd  u^o  (AT),  and  the  first-order  terms  are  ( l/M) 
Uci  (.V)  and  (1/M)  i/^,  (X). 

To  zeroth  order  in  the  core  C 1  the  velocity  i/ci.o  =  1 
throughout  the  region.  The  velocity  is  the  same  as  the  mov¬ 
ing  wall  throughout  this  core  region.  The  two  components  of 
the  current  density  to  zeroth  order  i.o  aid  j^ci.o  are  both 
equal  to  zero.  Also  the  potential  to  zeroth  order  throughout 
the  core  Cl  is  t^c  i.o  =  ^  throughout  the  core  area.  To 

zeroth  order  in  core  C  2  the  current  density  follows  the  exter¬ 
nal  magnetic  field  B„  and  therefore  y^c2.o  =  0.  and 
JriC2.o  =  —  cos  0((^o -f  ^)  [  see  Eq.  (17a)].  The  zeroth  order 
potentials  and  ui-q,^)  are  functions  of  q  and  ^  [see 

Eqs.  (7b)  and  (7c)  ] .  It  should  be  noted  that  to  zeroth  order 
there  is  a  jump  on  i/c  2  across  each  Hartmann  layer  at  the  top 
and  bottom  wall  [ see  Eq.  ( 1 7e)  j .  In  core  C  3  to  zeroth  order 
“c3.o  zero  throughout  the  subregion.  Thus,  the  core  C3 
flow  has  no  velocity  throughout  the  subregion.  Both  compo¬ 
nents  of  the  current  density  in  C  3,  y  ,c  3,0  and  y  ^  c  3.0  are  zero 
as  in  core  subregion  Cl .  To  zeroth  order  the  potential  ^c3,o  *s 
also  zero  throughout  the  region.  It  was  not  necessary  at  this 
point  in  the  work  to  obtain  the  Hartmann  layer  solutions 
adjacent  to  each  core  region  to  solve  for  the  core  region  flow 
profiles  to  zeroth  order. 

The  fundamental  mathematical  theory  was  developed 
for  the  free-shear  layers  /land  /2in  the  rectangular  chan¬ 
nel.  The  free-shear  layers  in  regions  f\  and  fl  in  coordi¬ 
nates  (x:,  t,  ^)  are  shown  in  Fig.  4.  The  corner  region  of  the 
free-shear  layer  / 1  between  the  insulating  wall  and  the  per¬ 
fect  conductor  is  at  (  /  =  0,  !■=  —  a  cos  6).  The  free-shear 
layer  fl  has  the  corner  region  between  the  insulating  wall 
and  the  perfect  conductor  at  (t  =  sec0,  ^  =  acos 
0  —  sin0).  In  this  work  we  developed  the  fundamental 
mathematical  theory  for  shear  layer  / 1 :  the  theory  for  shear 
layer  fl  can  be  performed  in  an  analogous  manner.  Alty'^ 
previously  studied  this  problem  of  the  free-shear  layer  in  a 
rectangular  channel  with  all  stationary  walls,  but  his  free- 
shear  layer  solution  is  wrong. 

The  basic  fourth  order,  partial  differential  equation  for 
the  potential  to  zeroth  order  in  / 1,  [Eq.  (47a)  ]  in 

the  free-shear  layer  was  derived.  The  solution  to  the  partial 
differential  equation  fits  five  different  boundary  conditions 
[see  Eqs.  (47b)-(47e)  j.  The  partial  differential  equation 
was  solved  in  this  work  by  using  the  mathematical  technique 
of  the  complex  Fourier  transform.  The  mathematical  solu- 
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tion  to  the  problem  using  this  technique  evolves  into  an  inte¬ 
gral  equation  which  can  be  solved  numerically  for  the  ze- 
roth-order  velocity  profile  in  shear  layer  /  I.  The  numerical 
solution  of  the  integral  equation  will  be  presented  in  a  future 
paper.  From  the  fundamental  theory  for  the  free-shear  layer 
the  basic  characteristics  of  the  free-shear  layer  were  deter¬ 
mined  and  discussed  in  detail.  There  has  been  experimental 
verification  of  a  similar,  but  not  exactly  the  same,  type  of 
free-shear  layer  as  presented  herein,  reported  by  Hunt  and 
Malcolm. ' ' 

The  basic  theory  was  derived  to  calculate  to  zeroth  or¬ 
der  the  total  current  per  unit  length  of  the  channel  in  core  C  2 
due  to  an  external  potential  applied  across  the  electrodes 
( i.e.,  perfect  conducting  walls ) .  Also,  the  expressions  for  the 
power  losses  per  unit  length  of  the  channel  were  derived  to 
zeroth  order  for  core  C  2.  The  derived  expressions  are  depen¬ 
dent  on  channel  aspect  ratio  a,  external  magnetic  field  angle 
6,  and  Hartmann  number  M.  No  current  is  carried  in  cores 
C 1  and  C  3  to  zeroth  order.  The  total  current  per  unit  length 
terms  neglected  in  the  two  free-shear  layers  and  the  two 
Hartmann  layers  along  the  insulating  walls  are  of  0(  Af  ~ ) 
and  0{M  ~ '),  respectively. 

In  this  paper  Figs.  9-14  show  linear  plots  of  the  mean 
voltage  at  the  center  of  the  perfectly  conducting  moving 
wall  versus  total  current  per  unit  length  of  the  channel  (// 
Aa).  The  potential  <f>„  is  referenced  to  the  perfect  conduct¬ 
ing  wall  at  the  bottom  of  the  channel,  which  is  assumed  to 
have  a  potential  of  <^  =  0.  The  linear  plots  show  that  for  the 
range  of  variables  studied  for  positive  mean  voltage  the 
total  current  per  unit  length  is  negative.  The  total  current  is 
thus  flowing  from  the  moving  wall  to  the  bottom  perfect 
conductor  in  the  duct  in  core  C2.  Similar  linear  plots  show 
that  for  negative  mean  voltages  <b„  the  total  current  is  posi¬ 
tive  and  is  flowing  from  the  bottom  stationary  conducting 
wall  to  the  top  moving  conducting  wall  in  core  C  2. 

The  following  general  conclusions  can  be  drawn  from 
this  range  of  calculated  data.  The  plots  show  that  as  the 
positive  mean  voltage  cb,„  increases,  as  the  aspect  ratio  of 
channel  a  increases,  and  as  the  external  magnetic  field  angle 
6  decreases,the  total  negative  current  per  unit  length  in¬ 
creases.  The  total  current  is  flowing  from  the  top  moving 
electrode  to  the  bottom  moving  electrode.  The  plots  also 
show  that  as  the  negative  mean  voltage  increases,  as  the 
aspect  ratio  of  the  duct  increases,  and  as  the  external  mag¬ 
netic  field  angle  6  decreases,  the  total  positive  current  per 
unit  length  increases.  The  total  current  is  flowing  from  the 
bottom  to  the  top  electrode.  It  should  be  noted  at  this  point 
that  positive  mean  potential  d)„  gives  a  larger  absolute  value 
total  current  per  unit  length  value,  |  (//Ax)  |,  than  the  corre¬ 
sponding  negative  potentials  in  a  similar  channel. 

Figure  1 5  shows  a  plot  of  nondimensional  resistance  per 
unit  length  of  channel  versus  angle  of  the  external,  homogen¬ 
eous  magnetic  field  B„.  For  this  data  the  conclusion  is 
reached  that  the.  resistance  of  the  channel  increases  with  in¬ 
creasing  external  magnetic  field  6  and  decreasing  channel 
aspect  ratio  a,  as  would  be  expected. 

Nondimensional  Joulean  power  losses  per  unit  length  of 
channel  versus  mean  potential  for  different  external 
magnetic  field  angles  6  and  aspect  ratios  a  are  plotted  in  Figs. 
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16  and  17.  For  limited  data,  power  losses  increase  with  in¬ 
creasing  positive  potential  or  decreasing  negative  potential 

and  with  decreasing  external  magnetic  field  angle  6. 

For  fully  developed  flow,  we  assume  that  there  is  only 
one  component  of  velocity  ux  and  that  all  variables  are  inde¬ 
pendent  of  X  and  time  t  *.  Therefore,  the  inertial  term  in  the 
Navier-Stokes  Eq.  ( la), 

is  identically  equal  to  zero  for  any  value  of  the  interaction 
parameter  N.  Fully  developed  flow  is  always  a  mathemat¬ 
ically  correct  solution  of  the  governing  equations,  but  it  may 
not  actually  occur  in  the  real  duct  flow.  Under  certain  cir¬ 
cumstances,  the  fully  developed  flow  in  unstable  with  re¬ 
spect  to  a  small,  three-dimensional,  time  varying  distur¬ 
bance.  The  disturbance  grows  until  a  different  stable  flow  is 
realized.  In  the  actual  flow,  all  three  components  of  velocity 
are  nonzero  and  the  flow  variables  depend  on  both  x  and  t  *. 
This  flow  may  involve  a  coherent  structure,  such  as  a  series 
of  vortices,  or  a  chaotic  structure,  such  as  turbulence.  In 
magnetohydrodynamic  flows,  the  magnetic  field  strongly  af¬ 
fects  the  instability  and  the  structure  of  the  periodic  flow. 

For  the  present  problem,  the  largest  velocities  occur  in 
the  free-shear  layers  / 1  and  fl,  so  that  inertial  effects  will 
first  become  significant  in  these  layers.  The  key  inertial  pa¬ 
rameter  for  these  layers  is  a  M  Ifa<  1,  then  the 

electromagnetic  body  force  overwhelms  inertial  effects,  the 
inertial  terms  in  Eq.  (la)  are  negligible,  even  in  the  high 
velocity  free-shear  layers,  and  fully  developed  flow  is  real¬ 
ized.  This  corresponds  to  the  statement  that  low  Reynolds 
number  flows  are  stable  in  ordinary  hydrodynamics.  We  can 
increase  a  by  increasing  f/„  or  by  decreasing  5„.  At  some 
critical  value  of  a,  fully  developed  flow  becomes  unstable 
and  there  is  a  transition  to  a  periodic  flow.  For  a  slightly 
above  a„,  there  is  a  velocity  fluctuation  wave  in  the  free- 
shear  layer.  This  wave  propagates  with  a  definite  wavelength 
and  wave  speed  around  the  circumference  of  the  current  col¬ 
lector  which  is  modeled  here.  At  first,  this  wave  would  repre¬ 
sent  a  propagating  undulation  in  the  free-shear  layer.  As  a  is 
increased  further,  the  undulations  grow  and  become  a  series 
of  vortices  whose  axes  are  parallel  to  B  in  the  plane  of  the 
free-shear  layer.  The  vortices  are  convected  around  the  cir¬ 
cumference  of  the  current  collector.  The  vortices  provide 
significant  momentum  transport  in  the^  direction.  This  mo¬ 
mentum  transport  causes  the  free-shear  layer  to  become 
thicker  and  its  maximum  mean  velocity  to  become  smaller. 
The  layer  thickness  grows  from  while  the 

mean  u  decreases  from  M''-  to  N''^.  However,  inertial  ef¬ 
fects  are  still  negligible  in  the  core  regions,  which  involve 
much  smaller  velocities,  so  that  the  present  core  solutions 
still  apply.  The  core  velocities  and  potentials  are  the  same,  so 
that  the  jump  in  potential  across  each  free-shear  layer  is  the 
same.  The  relationship  Uj  =  dd)^  /d^  for  each  free-shear  lay¬ 
er  comes  from  ohms  law,  so  it  is  unchanged.  The  new  period¬ 
ic,  inertial  free-shear  layers  with  0{N  ^ ''’)  thickness  must 
have  the  same  total  flow  at  each  t  =  constant  level  as  the 
original  thickness  layers. 

In  the  free-shear  layers  for  fully  developed  flow,  there  is 
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a  large  viscous  dissipation  because  the  velocity  gradient  dUf/ 
=  0{M).  After  transition  to  the  periodic,  inertial  free- 
shear  layer  structure,  the  viscous  dissipation  in  the  mean 
velocity  is  much  smaller  since 


dUf 

¥ 


=  0(N-'^)4:0(M). 


However,  the  mean  flow  also  loses  energy  to  the  vortices 
through  the  Reynolds  stress.  The  vortices  receive  energy 
from  the  mean  flow  and  lose  it  to  viscous  dissipation  and  to 
electromagnetic  damping  (Joulean  dissipation).  An  infi¬ 
nitely  long  vortex  with  rotation  around  a  straight  magnetic 
field  line  represents  a  generator  on  open  circuit:  The  radial 
component  of  vXB  is  exactly  cancelled  by  the  radial  electric 
field,  so  that  there  is  no  current  and  no  electromagnetic 
damping.  However,  in  the  present  problem,  each  magnetic 
field  line  intersects  one  or  two  perfect  conductors.  These 
perfect  conductors  provide  a  short  circuit  for  the  vortex  gen¬ 
erator.  Electric  current  circulates  through  the  conductor 
and  vortex  and  strongly  damps  the  vortex.  If  a  is  increased 
further,  nonlinear  vortex  interactions  lead  to  a  wide  spec¬ 
trum  of  vortex  sizes,  i.e.,  to  turbulence.  However,  the  turbu¬ 
lence  is  two-dimensional  since  any  eddy  whose  axis  is  not 
parallel  to  B  is  strongly  damped."  The  value  of  a  when  each 
transition  occurs  depends  critically  on  the  geometry  and 
conductivity  of  the  duct  walls. 
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APPENDIX  A:  THICKNESS  OF  LAYERS  IN  FULLY 
DEVELOPED  MAGNETOHYDRODYNAMIC  CHANNEL 
FLOWS 


The  discussion  herein  presents  a  fundamental  general 
magnetohydrodynamic  mathematical  treatment  for  deter¬ 
mining  the  order  of  magnitude  of  boundary  layer  thickness 
in  fully  developed  laminar,  magnetohydrodynamic  channel 
flows. 

Equations  for  fully  developed  magnetohydrodynamic 
channel  flows  with  a  uniform  magnetic  field  are  as  follows: 


Jv  =  - 


\dv^ 
d<f) 
dr}  ’ 


+ 


d^u\ 

di^j’ 


(Ala) 

(Alb) 


d<t>  , 

^-1-^  =  0. 

dv  di 


(Ale) 


(Aid) 


Here  (x,  r),  are  Cartesian  coordinates  with  the  velocity  « 
in  the  x  direction  and  the  magnetic  field  in  the  rj  direction. 
Introduce  Eq.  (Ale)  for  the  u  in  Eq.  (Ala)  to  obtain 


-'  3  ^2  +  gtz)- 


.dT]-di  ■  di^  '  dt}-  ’  di 
Now  we  differentiate  Eq.  ( A2)  with  respect  to 

-2  f  d*4>  j  d*^  ^ 


dr}-di  di 


dl^J' 


di  Kdrj^df-  d^* 

Substitute  Eq.  (Alb)  intoEq.  (Aid)  to  obtain 

'd^  dr}  dr}-  ■ 

Substitute  Eq.  (A4)  for  each  into  Eq.  (A3) 


(A2) 


(A3) 


(A4) 


dr} 

d-4> 

dr}- 


drj^df- 


d^*  drf  dr}-d^ 


(A5) 


Now  let  us  consider  a  wall  in  which  the  unit  normal  y  is 
at  an  angle  0  to  the  magnetic  field  direction  (see  Fig.  18). 
The  unit  vector  -ii,  the  first  derivative  d /dr},  and  the  second 
derivative  can  be  written  as  in  terms  of  Cartesian  coordinates 


•q  =  cos  0y  -f  sin  0i, 

—  =  cos  6 - h  sin  p  — , 

dr}  dy  dz 


(A6a) 

(A6b) 
az 

=  cos’  6  -4-  +  2  sin  0  cos  6  — - h  sin’  6  . 

dr}-  dy-  dydz 

(A6c) 

Therefore  in  ix,y,  z)  coordinates,  Eq.  (A5)  is  expressed  as 

2  21^V  ,  1  ■  21  21  ,  •  221^V 

cos  -I-  2  sm  0COS  6-r^  +  sivr  0-  — 


dydz 


dr 


=  "-’(4^  +  2 

Vol/ 


dfd^ 

^  dz*) 

(A7) 


Here  it  noted  that  V"*  is  unchanged  by  rotating  the  axis. 
Equation  ( A7 )  is  a  fourth  order  elliptic  equation  so  that  two 
boundary  conditions  must  be  described  at  each  boundary. 
For  a  magnetohydrodynamic  problem,  the  velocity  is  usual¬ 
ly  u  =  0  and  either  the  potential  (f>  is  given  at  a  perfect  con¬ 
ductor  or  jy  could  be  given  at  a  finite  conductor  as 
jy  =  C(d^<p/d^)  where  C  =  {a^t/crL)  is  the  wall  conduc¬ 
tance  ratio,  where  cr^  and  t  are  the  conductivity  and  thick¬ 
ness  of  the  wall,  respectively.  Therefore,  let  us  assume  that 
there  are  two  boundary  conditions  at  the  wall  at>>  =  0.  Now 
consider  the  asymptotic  solution  as  M-* « .  Unless  some  de¬ 
rivative  becomes  very  large  <f>  is  governed  by  the  “reduced 
equation”: 


cos’ 0 -f  2  sin  0  cos  0 -1- sin’ 0 4-^ .  (A8) 

dy^  dydz  d^ 

The  reduced  equation  is  only  second  order,  so  it  cannot  satis¬ 
fy  both  boundary  conditions  aty  =  0.  Therefore,  there  must 
be  a  boundary  layer  adjacent  to  the  wall  where  the  normal 
derivative  d  /dy  is  sufficiently  large  that  at  least  one  of  the 
fourth-order  derivatives  on  the  right-hand  side  of  Eq.  (A7) 
is  comparable  to  the  second-order  derivatives  on  the  left- 
hand  side.  Then,  near  the  wall,  the  governing  equation  again 
becomes  fourth  order  and  can  satisfy  both  boundary  condi¬ 
tions  at;’  =  0.  The  core  region  is  obviously  governed  by  the 
“reduced  equation.” 
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We  now  change  y  to  stretched  boundary  layer  coordi¬ 
nates. 


APPENDIX  B:  MATCHING  PRINCIPLE  BETWEEN 
HARTMANN  LAYER  AND  CORE  REGION 


r=5"'y,  (A9a) 

y=(5)T.  (A9b) 


i.  =  .^A  =  5->  A 

dy  dy  dY  dY 


(A9c) 


Equation  ( A7)  is  expressed  in  stretched  boundary  layer  co¬ 
ordinates  as 


cos=  d{8--) 


+  2  sin  0  cos  0(5  ') 
dY~ 


d-i 

dYdz 


-t-  sm'  6  — 
dr 

=  A/-^((5-')|^-F2(5-’) 


0V 

dY-dr 

dz*  } 

(AlO) 


The  first  term  on  the  right-hand  side  of  Eq.  (AlO)  is  always 
much  larger  than  the  other  two  terms.  There  are  three  cases 
for  various  5.  In  case  I,  cos0  is  not  small  [i.e.,  O<0<0„ 
<  (7r/2)  ],  so  the  magnetic  field  has  a  finite  component  nor¬ 
mal  to  the  wall.  Then  5  =  Af  “ '  and  the  basic  governing 
equation  is 


cos-  0 


dY- 


d^ 

dY*  ■ 


(All) 


This  is  the  governing  equation  for  all  Hartmann  layers  for 
walls  where  there  is  a  finite  normal  component  of  B.  The 
basic  solutions  are 


A,(z),  A2(z)Y,  ^,(z)exp(  -  ]cos0|y), 

A^iz)e\p(  4-  |cos0|y), 

where  the  i4,  (z)  are  arbitrary  functions.  Thus,  all  Hartmann 
layer  solutions  are  composed  of  linear  combinations  of  these 
four  solutions. 

In  case  II  the  angle  6  =  rr/1  and  the  magnetic  field  is 
exactly  parallel  to  the  wall.  Then  5  =  M  ~  and  the  gov¬ 
erning  equation  is  of  the  mathematical  form 


The  matching  principle  was  originally  stated  by  Lager- 
strom-"  -'  in  the  early  1560s,  which  was  proven  to  be  rigor¬ 
ous  by  O’Malley^-  -  '  in  the  late  1960s.  The  principle  states 
that 

n  inner(/n  outer)  =  m  outer(n  inner)  (B1 ) 

where  outer  =  core,  and  inner  =  boundary  layer.  We  will 
illustrate  with  an  example  2-inner  (2-outer)  =  2-outer  (2- 
inner)  where  (n  =  m  =  2). 

The  outer  (core)  expansion  is 

Uc(x\M)  =Ucoix)  +M~'uciix)  +M~-Uc2ix)  +  •••, 
2-outer  =  truncated  after  first  two  terms.  ( B2a) 

Thus, 

2-outer  =  Uco{-*^) -b  (B2b) 

Now,  inner  (n  outer)  =  n  outer  rewritten  in  inner 
(boundary  layer)  variables  so  that 

inner  (2-outer)  =  Uco  +  M~  'ua  {M  “  'A'). 

(B3) 

This  is  not  a  legitimate  asymptotic  expansion  because 
the  coefficient  functions  Uc,{M~*X)  depend  on  A/.  We  must 
change  this  to  asymptotic  expansion  before  we  go  on. 

For  AT  fixed,  AT  “'AT—O  as  A/— 00.  (B4) 

Therefore,  we  use  expansions  for  small  argument,  i.e.,  Tay¬ 
lor  series: 


Uco iM ■  'X)  =  Uco  m+M- 

dx 


1 

+  1a/-’a'’ 
6 


d'^UcoiO) 

dx^ 

d^UcoiO) 

dx^ 


(B5) 


Also  expanding  {M  'A()  in  a  Taylor  series  results  in  the 
following  expansion  for  the  truncated  series 


d^  _ 

dz^  dY*  ■ 


(AI2) 


This  is  the  governing  equation  for  “side  layers.”  These  are 
boundary  layers  on  walls  parallel  to  the  magnetic  field,  and  it 
is  the  governing  equation  for  free-shear  (or  interior)  layers 
that  lie  along  the  magnetic  field,  such  as  the  free-shear  layers 
being  studied  within  this  paper. 

Case  III  deals  with  the  problem  when  cos  0  =  aM 
and  8  =  where  a  is  an  0(  1 )  parameter  in  the  prob¬ 

lem.  Then  the  governing  Eq.  (AlO)  becomes 


inner  (2-outer) 


=  “co  (0)  +  M  ' 


-f  Af- 


t/Uci  (0) 
dX 


rfMco(O) ^ 
(O)-fAr— 


(B6) 


This  is  now  a  valid  asymptotic  expansion  in  the  form 

2MX)M-'.  (B7) 

I 

Therefore, 


■,  d~(/>  .  d~4>  ,0V  0V 

a  — —  -t-  la  — —  H - ^  ^  . 

dY-  dYdz  dz-  dY* 


(A13) 


This  is  the  equation  for  a  wall  that  is  nearly,  but  not  exactly, 
parallel  to  the  magnetic  field.  The  governing  equation  for 
case  III  becomes  that  for  case  II  as  a— 0  and  becomes  that 
for  case  I  as  a  —  oo ,  if  we  redefine  5. 


2-inner  (2-outer)  =  truncated  after  second  term, 

2-inner  (2-outer) 

=  Uco (0)  +  A/  - '  (nci  (0)  -b  X  )  •  (B8a) 

The  inner  expansion  can  be  expressed  as 
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inner  =  ttfco  W  +  3/ ~ 'a* ,  (A')  +  A/ ~-Ut2  W, 

2-inner  =  a*o(Ar)  +  3/“'ai,  (X), 

outer  (2-inner)  =u^o(Mx)  +3/“'a*,  (A/x).  (B8b) 

This  is  not  a  legal  expansion.  Thus,  we  need  to  know 
how  u^,(X)  behaves  as  Af—  oo . 

In  general,  we  have  solution  forms  from  Hartmann 
boundary  layer  equations,  so  we  know  how  the  solution  be¬ 
haves.  Assume,  for  some  a  >  0 

a*„(A')  —Aq  -f  /4..  exp(  —  aX)  as  Af—  oo  (B9) 

and 

U/,,  (x) -‘A,  +  B,X  +  e\p(  —  aX)  asAf—oo. 

By  inspection  the  following  relationships  can  be  ob¬ 
tained  for  the  Hartmann  layer  equations: 

lim  (a^o)  =^().  (BlOa) 

X  ■  K 

lim  (a^i  —  5|A')  = /!,.  (BlOc) 

A  • « 

Note  as  Af  —  00,  exp  (  —  aX)  4X  ~  "  for  any  «.  The  expan¬ 
sion  in  terms  of  these  expressions  is 

outer  (2-inner) 

=  A„  +  A^  exp(  —aAfx)  +Af~‘[A,  +B,Afx 
-i- B,.  exp(  —  aAfx)].  (Bll) 

Sinceexp(  —aMx)  =exp(  — aA')— OasA'— oo  2outer(2- 
inner)  =  A„  -f  +  M"' A,  +  asymptotically  exponen¬ 
tially  small  terms. 

Using  the  matching  principle  (Eq.  (Bl)]  as  Af— oo  in 
the  boundary  layer  we  obtain 

A„  +  B,M  -'X  +M  -  'A, 

=  Uc„(0)  +  Af-'uc,(0)+M 

ax 

(B12) 

Therefore,  the  zeroth  and  first  order  matching  condi¬ 
tions  between  the  core  region  and  Hartmann  layer  are 

/t„  =  «co(0)  =  lim  [«M)(-^)].  (B13a) 


/I,  =  w,.,(0)  =  lim  [a^iCAT) -fl.Ar],  (B13b) 

X  •  oc 


duci)(0) 

dx 


lim 

A  « 


dui,t(X) 

dX 


(B13c) 


APPENDIX  C:  MATCHING  PRINCIPLE  BETWEEN  THE 
FREE-SHEAR  LAYER  AND  HARTMANN  LAYER 

The  free-shear  layer  is  separated  from  each  wall  by  a 
Hartmann  layer  with  0(M~')  thickness.  Since  the  free- 
shear  layer  is  <?(  A/  thickness  and  .3/  ~*4M  these 

Hartmann  layers  appear  to  the  free-shear  layer  to  have  zero 
thickness.  That  is,  the  Tree-shear  layer  solution  applies  to 


within  an  arbitrarily  small  distance  from  at  wall  and  then 
there  is  a  jump  across  a  zero  thickness  layer  to  satisfy  the 
wall  conditions.  When  we  stretch  the  normal  coordinate  to 
the  wall  by  M,  we  can  look  at  the  structure  of  the  Jump  across 
the  “zero-thickness”  layer  which  has  been  stretched  to  have 
0(1)  thickness.  However,  since  the  Hartmann 

layer  sees  the  free-shear  layer  to  be  at  infinity.  Specifically, 
the  value  which  the  free-shear  layer  “sees”  at  an  arbitrarily 
small  distance  from  the  wall  appears  to  the  Hartmann  layer’s 
stretched  coordinate  to  be  infinitely  far  from  the  wall. 

The  free-shear  layer  “sees”  horizontal  boundaries  at 
r  =  0  and  t  =  sec  6,  even  though  walls  are  not  horizontal 
here. 


The  coordinate  transformations  are  (see  Figs.  2  and  4) 

U  +  a)  cos  9  +t  sin  9, 

(Cla) 

y  =  tcos9  —  ^sm  9, 

(Clb) 

t=y  cos  9  +  {z  +  a)sin  9, 

(CIO 

{z  +  a)cos  —  j'  sin  9. 

(Cld) 

Let  us  now  consider  the  matching  between  the  free- 

shear  layer  and  insulated  Hartmann  layer.  The  stretched  co¬ 
ordinates  are 

Y=M'>^y,  y  =  M-''^Y  (C2a) 

Z  =  M{z  +  a),  z=  —  a  +  M~'Z.  (C2b) 

The  Hartmann  layer  solution  can  be  expressed  as 

U/f  —  M  *^~Uho(Y,Z) U/i )  (Y,Z).  (C3) 

The  free-shear  layer  coordinate  and  stretched  coordinate 
are,  respectively, 

/  =  /,  (C4) 

The  velocity  in  the  free-shear  layer  can  be  expressed  as 
«/,  =  2-outer  =  (C5) 

For  inner  (2-outer),  we  need  t  and  f  in  terms  of  Y  and  Z. 
Note  that  t  can  be  expressed  in  terms  of y  and  z  as 

/  =>'cos  0 -i- (z -I- a)sin  0.  (C6) 

Substituting  Eqs.  (C2a)  and  (C2b)  for  stretched  coordi¬ 
nates  results  in  the  expression 

t  =  cos0(A/-‘'-T) -f  sin0(A/-'Z).  (C7a) 

Also,  ^  is  expressed  as 

f  =  A/''"|-  =  A/''-(A/-'Zcos0-A/-''=ysin  6) 

=  (cos(9A/-'Z- Tsin^).  (C7b) 

Now  the  free-shear  layer  velocity  can  be  expressed  as 

inner  (2-outer) 

=  M  I  (cos  eM  -'^-Y  +  sin  0M  '  'Z). 

(cosA/  ''’Z-ysin0)]-fMy, ,((  ),(  )]-)-•••. 

(C8) 

The  terms  inner  and  outer  are  defined  in  Appendix  B.  Next 
we  need  to  use  a  double  Taylor  series  expansion: 
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where 


«yj  (r,f)  =  «<;{  (0.  -  ysin  d)  +  (0,  -  rsin  O) 

at 

+  (f  +  rsin  0)  — ^  (0,  -  ysin  0),  (C9a) 

_ I 


(f  +  ysin  0)=co&0M-''^Z. 

Substituting  in  Eq.  (C8)  results  in  the  expression 


2-inner(2-outer)  —  ysin 5)  +  ico&0M~^'^Y  +  %m0M~'Z) — ^  (0,  —  ysin0) 

\  dt 


Keeping  only  M  '^^and  0(  1 )  terms  for  2-inner  (2-outer  expression; 


duf. 


2-inner  (2-outer)  =  —  ysin  0)  -b  cos  0Y — —  (0,  —  ysin  0) 

dt 

^„(0) 

+  cos  0Z  — —  (0,  —  y  sin  0)  -f  «/,’(0,  —  y  sin  0). 

When  we  do  the  matching  all  conditions  on  «yi  will  be  at  /  =  0,  so  it  sees  a  horizontal  wall 


(C9b) 


^y(O)  \ 

-t-  (cos  0M  - '/^Z)  — ^  (0,  -  ysin  6>)  -t-  •  •  •  +  «y,’(0,  -  ysin  0)  4-  (cos  0M -''^Y 

J 

-f  sin  0M - 'Z)  — ^  (0,  -  ysin  0)  +  cos  0M “ '^^Z — ^  (0,  -  ysin  0)  +  ---. 
dt  d^ 


(CIO) 


(Cll) 


1 -inner (1 -outer)  =  1 -outer ( 1 -inner):  lim  Ufjo(Y,Z)  =  «‘n’(0,  —  ysin^). 

Z-co 


(C12) 
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In  a  previous  paper  the  authors  initiated  studies  of  fully  developed  laminar  liquid-metal  flows, 
currents,  and  power  losses  in  a  rectangular  channel  with  a  moving  perfectly  conducting  wall 
and  with  a  skewed  homogeneous  external  magnetic  field  for  high  Hartmann  numbers,  high 
interaction  parameters,  low  magnetic  Reynolds  numbers,  and  different  aspect  ratios.  The 
channel  had  insulating  side  walls  that  'vere  skewed  to  the  external  magnetic  field,  while  the 
perfectly  conducting  moving  top  wall  with  an  external  potential  and  the  stationary  perfectly 
conducting  bottom  wall  at  zero  potential  acted  as  electrodes.  These  electrodes  were  also 
skewed  to  the  external  magnetic  field.  A  mathematical  solution  was  obtained  for  high 
Hartmann  numbers  by  dividing  the  flow  into  three  core  regions,  two  free  shear  layers,  and  six 
Hartmann  layers  along  the  channel  walls.  The  free  shear  layers  were  treated  rigorously  and  in 
detail  with  fundamental  magnetohydrodynamic  theory.  The  previous  work,  however,  left  the 
solution  for  the  velocity  profiles  in  terms  of  a  complex  integral  equation  which  was  not  solved. 

In  the  present  work  the  integral  equation  is  solved  numerically  by  the  method  of  quadratures 
to  give  the  velocity  profiles,  viscous  dissipation  and  Joulean  losses  in  the  free  shear  layers.  In 
addition,  expressions  for  the  viscous  dissipation  in  the  six  Hartmann  layers  are  presented.  The 
best  approximation  to  the  viscous  dissipation  in  the  channel  is  the  sum  of  the  0(Af 
contributions  from  the  two  free  shear  layers,  the  contributions  from  the  two 

Hartmann  layers  separating  the  free  shear  layers  from  the  insulators,  and  the  0{M) 
contributions  from  three  of  the  Hartmann  layers  separating  core  regions  from  the  walls.  The 
best  approximation  to  the  Joulean  power  losses  in  the  channel  is  the  sum  of  the  0(3/^) 
contribution  from  the  central  core  region  which  carries  an  0(  1 )  current  between  the  electrodes 
and  the  contributions  from  the  free  shear  layers.  The  expressions  for  the  viscous 

dissipation  and  Joulean  losses  in  each  region  involve  the  products  of  universal  constants, 
electrical  potentials  and  geometric  factors.  The  theoretical  magnetohydrodynamic  model 
presented  here  was  developed  to  provide  data  to  help  in  the  design  of  liquid-metal  current 
collectors. 


I.  INTRODUCTION 

Advanced  homopolar  electrical  machinery  is  being  de¬ 
veloped  for  many  new  applications.  These  advanced  ho¬ 
mopolar  machines  require  current  collectors  for  transport¬ 
ing  high  current  at  low  voltages  between  rotating  and 
stationary  members  of  the  machines  with  minimal  losses  and 
maximal  operational  stability.  Therefore,  liquid  metals  are 
often  used  in  the  sliding  electric  contact  region,  instead  of 
more  conventional  brush  technology.  Generally,  these  ma¬ 
chines  have  large  magnetic  fields  in  the  region  of  the  current 
collector,  which  produce  Lorentz  forces  on  the  liquid  metal 
which  can  increase  power  losses  and  lead  to  flow  instabili¬ 
ties.  Thus,  to  obtain  current  collector  design  parameters,  one 
is  interested  in  liquid-metal  channel  flow  problems  with  ap¬ 
plied  external  magnetic  fields  and  with  boundary  conditions 
containing  combinations  of  moving  and  fixed,  conducting, 
or  insulating  walls.  These  boundary  conditions  then  approx¬ 
imately  correspond  to  a  rotor,  stator,  or  sidewall  of  the  cur¬ 
rent  collector  and  the  liquid  metal  corresponds  to  the  sliding 
electric  contact. 


In  contrast  to  the  rather  extensive  magnetohydrodyna¬ 
mic  (MHD)  duct  flow  literature,’”’  investigations  treating 
moving  conducting  walls,  such  as  discussed  by  Hughes  and 
Young,”  are  rare.  Brown,  Reilly,  and  Sondergaard"  have 
studied  fully  developed,  viscous  liquid-metal  flows  and  pow¬ 
er  losses  in  a  rectangular  channel  with  a  moving  conducting 
wall  in  a  uniform,  external  magnetic  field  for  moderate  Hart¬ 
mann  numbers.  The  rectangular  channel  had  insulating  side 
walls  parallel  to  the  magnetic  field  and  a  perfectly  conduct¬ 
ing  moving  top  wall  and  stationary  bottom  wall  perpendicu¬ 
lar  to  the  field.  A  previous  paper"  by  the  authors  treated 
rigorously  and  comprehensively  this  same  problem  with  two 
major  differences.  The  homogeneous  transverse  magnetic 
field  was  changed  to  a  high  strength,  skewed  magnetic  field, 
and  an  external  potential  was  applied  to  the  conducting 
moving  wall,  thus  making  the  channel  transport  current. 
Skewed  magnetic  fields  are  present  at  the  collectors  in  some 
machine  designs  and  external  potentials  on  the  electrodes 
correspond  to  the  generated  or  motored  voltage  of  the  ma¬ 
chines.  These  changes  complicate  the  channel  problem  by 
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FIG  1  Rectangular  channel  with  skewed  external  magnetic  field  and  with 
moving  and  stationary  perfect  conductors 


removing  symmetry.  The  present  paper  should  be  consid¬ 
ered  as  part  two  of  the  previous  paper.*’ 

The  purpose  of  the  previous  paper*’  was  to  present  com¬ 
prehensive,  two-dimensional  analytical  solutions  for  the 
flow  velocity  and  electrical  potential  in  different  regions  of 
the  channel  (see  Fig.  1).  In  the  three  core  regions,  the  singu¬ 
lar  perturbation  problem  for  the  velocity  and  electric  poten¬ 
tial  was  solved  to  zero  order  using  matched  asymptotic  ex¬ 
pansions.  The  solutions  in  the  three  core  regions  enabled  the 
investigation  of  the  influence  on  the  velocity  and  electrical 
potential  of  the  insulating  sidewalls  and  an  external  skewed 
high  magnetic  field.  Furthermore,  analytic  expressions  were 
derived  for  the  quantities  of  technical  interest  in  the  cores  to 
zero  order.  These  quantities  were  the  load  currents  between 
the  electrodes  (i.e.,  moving  and  stationary  conducting 
walls),  Joulean  and  viscous  dissipation,  and  channel  resis¬ 
tance  with  skew  angle.  The  parameters  of  variation  were 
applied  skewed  external  magnetic  field,  fluid  electrical  c  n- 
ductivity  and  viscosity,  and  channel  aspect  ratio. 

The  fundamental  magnetohydrodynamic  theory  for  the 
two  free  shear  layers  was  also  derived  to  zero  order.  How¬ 
ever.  the  cxpre.ssions  for  the  velocity  and  electrical  potential 
were  left  in  a  form  requiring  the  solution  of  an  integral  equa¬ 
tion. 

The  work  begun  in  our  earlier  paper*'  is  combined  and 
completed  here.  We  have  calculated  the  velocity  profile,  vis¬ 
cous  dissipation,  and  Joulean  power  losses  in  the  free  shear 
layers  by  numerically  solving  the  integral  equation  with  the 
method  of  quadratures.  Furthermore,  analytic  expressions 
for  the  velocity  profiles  and  viscous  dissipation  in  the  six 
Hartmann  layers  m  the  channel  are  pre.sented. 

In  treating  the  magnetohydrodynamic  flow  due  to  elec¬ 
tric  current  between  two  circular  electrodes  embedded  in 
two  parallel  insulators  with  a  transverse  magnetic  field. 
Alty*'  and  Hunt  and  Malcolm*'*  ''  treated  a  similar  free 
shear  layer,  but  their  approach  is  quite  different  from  ours. 

II.  NONOIMENSIONAL  MAGNETOHYDRODYNAMIC 
EQUATIONS 

The  nondimensional  magnetohydrodynamic  equations 
for  laminar  duct  flow  in  an  external,  homogeneous  magnetic 


field  at  an  angle  6  with  the  vertical  can  be  written  as 


)V-'(vV)v=  -V/*H-jXb-f  (la) 

j=  _V(^-t-vXb,  (lb) 

V-v  =  0,  (Ic) 

V.j  =  0,  (Id) 

where 


N  =  interaction  parameter,  the  ratio  of  ponderomotive 
force  to  inertial  force  =  aBlL  /pU„;  (2a) 

M  —  Hartmann  number,  the  positive  square  root  of 
the  ratio  of  the  ponderomotive  force  to  the  viscous 
force  =  I5„((r//r/)'^^  (2b) 

V  =  fluid  velocity  vector  normalized  by  U„  =  v*/C/„;  (2c) 
P=  pressure  normalized  by  aUifilL  =  P*/aU,^lL; 


(2d) 

j  =  electric  current  density  vector  normalized 

by  =  j*/CTC/,>5„;  (2e) 

=  electrical  potential  normalized  by 

=  (2f) 

V  =  gradient  normalized  by  L  =  LV*;  (2g) 

b  =  unit  vector  parallel  to  external  magnetic 

field  B„(B„  =  5„b);  (2h) 

(x,y,z)  =  unit  vectors  in  Cartesian  coordinates 
(.see  Fig.  1), 

b  =  cos  ^ -t- sin  ^z;  (2i) 


where  an  asterisk  denotes  a  dimensional  quantity.  Here  cr 
and  Pf  are  the  electrical  conductivity  and  viscosity  of  the 
liquid  metal,  while  L  is  the  distance  between  the  perfectly 
conducting  electrodes,  t/„  is  the  velocity  of  the  moving  per¬ 
fectly  conducting  top,  and  is  the  strength  of  the  external 
magnetic  field. 

It  is  assumed  during  this  work  that  the  magnetic  Reyn¬ 
olds  number  /?„,  =  U^^Lap  (i.e.,  ratio  of  induced  magnetic 
field  to  external  magnetic  field)  is  much  less  than  one,  where 
p  is  the  liquid’s  magnetic  permeability. 

In  the  Cartesian  system  of  dimensional  coordinates  (see 
Fig.  1 ),  the  (x.y.r)  coordinates  are  normalized  by  L.  The 
perfectly  conducting  sliding  wall  is  at  a  constant  velocity  of 
U„\  a\y=  1,  the  stationary  perfectly  conducting  wall  is  po¬ 
sitioned  at  y  =  0.  The  insulating  walls  are  at  z  =  ±a,  where 
2a  is  the  aspect  ratio  of  the  channel.  The  external  magnetic 
field  angle  0  is  represented  in  terms  of  the  field  components 
as 

(?  =  arctanffly^,. ),  (3a) 

and  the  magnitude  of  the  external  magnetic  field  in  terms  of 
the  field  components  is 

B„=(Bl  +  Biy‘-.  (3b) 


Equation  ( la)  is  the  dimensionless  Navier  Stokp  equa¬ 
tion  with  an  external  ponderomotive  body  force  jxb.  Equa¬ 
tion  (lb)  is  the  dimensionless  expression  for  the  electric  cur¬ 
rent  density  induced  in  the  channel.  Equation  (Ic)  is  the 
dimensionless  expression  for  the  fluid  incompressibility. 
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(7b) 


moving  porfoci  conductor 


FIG  2  Rectangular  channel  with  coordinate  systems  (^, >>,2)  and 


Equation  (Id)  represents  the  conservation  of  charge. 

A  new  rotated  nondimensional  coordinate  system 
(x,r;,g )  is  defined  with  the  j/  axis  parallel  to  the  external  field 
B„.  Therefore,  the  unit  vectors  1]  and  b  are  identical  (see  Fig. 
2),  and  x  is  the  unit  vector  in  the  axial  direction  along  the 
channel.  In  our  new  rotated,  right-handed  coordinate  sys¬ 
tem,  the  straight  lines  describing  the  top,  bottom,  left,  and 
right  walls  are  represented  by  the  following  equations: 


T]  =  sec  6  +  (tan 

for  -  a  cos  0  —  sin  cos  0  —  sin  0,  (4a) 

rj  =  (tan  6)^  for  —  a  cos  cos  0,  (4b) 

?7  =  —  a  CSC  <9  —  (cot  0}^ 

for  -  acosd  -  sin  —  a  cos  6,  (4c) 

7/  =  a  CSC  0  -  (cot  6)^ 

for  a  cos  0  -  sin  cos  6.  (4d) 

It  should  be  noted  that  the  unit  vectors  in  the  two  sys¬ 
tems  nave  the  following  relationships; 

z  =  sin  (^  +  cos  0 1,  (5a) 

y  =  cos  (Hj  -  sin  0 1,  (5b) 

In  our  rotated  coordinate  system  (x,r],^)  the  physically 
realizable  variables  for  fully  developed  flow"’  can  be  ex¬ 
pressed  as 

v  =  «(7/,s)x,  (6a) 

P=P(7].^)  (no  axial  pressure  gradient),  (6b) 

(j>  =  (6c) 

j  +7i  (6d) 


Substituting  the.se  variables  into  Eqs.  ( la)-(  Id)  results  in 
the  following  system  of  magnetohydrodynamic  equations  in 
terms  of  velocity  and  electrical  potential.  The  electrical  po¬ 
tential  and  flow  velocity  are  variables  that  can  be  measured 
experimentally. 


7.  =  4/ 


,  d-u\ 
\drf'^  df-} 


(7a) 


;  - 

dv  ’ 


Ja  = 


dr)  di 


(7c) 


and  F  is  a  constant. 

Here  we  assume  that  the  Hartmann  number  M  is  large. 
For  1 ,  the  interior  of  the  channel  is  subdivided  into  three 

core  regions  with  all  0(  1 )  derivatives,  two  free  shear  layers 
which  have  0(M~''^)  thicknesses  and  which  separate  adja¬ 
cent  core  regions,  and  Hartmann  layers  which  have 
0(M~')  thicknesses  and  which  separate  cores  or  free  shear 
layers  from  adjacent  walls.  In  each  core  region, 


4  =  J’l  =J’,  ^  =  -  viv  (i)> 


(8) 


neglecting  0(M  ~~)  terms.  The  jump  in  electrical  potential 
or  normal  current  density  across  a  Hartmann  layer  is 
0(M  '' ).  The  integration  functions (^)  and  i{i(^)  in  each 
core  are  determined  by  the  electric  potential  of  either  perfect 
conductor  or  by  the  boundary  condition  =  0  at  either  in¬ 
sulator,  neglecting  0(M  ')  terms.  This  completely  deter¬ 
mines  the  core  solutions,  which  we  presented  in  our  previous 
paper.’’  The  electrical  potentials  in  adjacent  cores  are  not 
equal  at  I"  =  —  a  cos  0  and  at  ^  =  a  cos  $  —  sin  8,  and  the 
free  shear  layers  at  these  planes  must  accommodate  these 
jumps  in  the  electrical  potential.  The  dimensionless  poten¬ 
tials  at  the  fixed  and  moving  conductors  are  zero  and 
(fl^d  +  i),  respectively. 


ill.  FREE  SHEAR  LAYER  ANALYSIS 

The  free  shear  layer  in  region /I  is  shown  in  Fig.  3.  The 
potential  at  the  left  is 

<i>c\  =<i>ii  —  a  cos  6,  (9a) 

and  the  potential  at  the  right  is 

=  cos  9((f>„  —  a  cos  6)(r}  a  sin  8).  (9b) 

The  width  of  the  layer  is  0(M  '  '^•)  and  the  top  and  bottom 
are  at  7/  =  sec  0  —  a  sin  8  and  7/  =  —  a  sin  0,  respectively. 


FIG,  3.  Free  shear  layer /I  separating  cores  cl  and  c2 
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We  shall  work  out  the  fundamental  mathematical  theory  for 
shear  layer/1;  the  theory  for  shear  layer/2  can  be  performed 
in  an  analogous  manner.  In  order  to  develop  the  theory  for 
/I,  the  stretched  coordinate  f  will  be  defined  as 

^  —  +  acos  6),  —  acos  0 

(10a) 


where  the  boundary  condition  from  matching  core  c2  is 

— 1)  as  ^—00,  (15b) 

the  boundary  condition  from  matching  core  cl  is 

«l>-.0  asf--oo,  (15c) 

the  boundary  condition  at  the  moving  perfect  conductor  is 


—  =  =  —  .  (10b) 

di  di  d; 

The  fundamental  magnetohydrodynamic  Eqs.  (7a)- 
(7c)  in  the  free  shear  layer  /I  for  fully  developed  flow  in 
stretched  coordinate  f  and  regular  coordinate  rj  are 


(11a) 


dt} 

4= 


(11b) 


%L  +  m''A  =  o_ 
dv 

The  first  terms  in  the  series  solutions  in  the  free  shear  layer/1 
are 


(12a) 

7^(77,^)  =A/-''V<n.o(’7.f)- 

(12b) 

u(77,p  =A/'^‘«/i.o(’7.^). 

neglecting  the 0{M  0(M  ~ ' ),  0{M  and  0(  1 ) 

terms,  respectively. 

From  the  magnetohydrodynamic  Eqs.  (I la)  and 
(lib),  the  following  relationships  can  be  developed: 


“/I.0  —  — 

dt] 

4'f\.o 

b 

II 

' 

11 

(13a) 


(13b) 


The  mathematical  solutions  to  the  differential  Eq.  ( 1 3a)  and 
the  boundary  values  of  Eqs.  (9a)  and  (9b)  are  assumed  to  be 
in  the  form 


<^/i.o  =<^o-acos0 +<l>(/,f),  (14a) 

where 


<1>  =  0  atr  =  sec0  for— oo<f<oo,  (15d) 

the  boundary  condition  at  the  fixed  perfect  conductor  is 

0  =  —  (<i„  —  a  cos  0)  at /  =  0  forO<f<oo, 

(15e) 

and  the  boundary  condition  at  the  Hartmann  layer  separat¬ 
ing  the  free  shear  layer  from  the  insulating  wall  is 

<1> 

,  +  —  =  0  at  /  =  0  for— oo<^<0.  (15f) 

The  boundary  conditions  at  the  perfect  conductors  apply 
because  the  jumps  in  electrical  potential  across  the  Hart¬ 
mann  layers  between  the  free  shear  layer  and  perfect  conduc¬ 
tors  are  0{M  “ ' ).  The  boundary  condition  Eq.  ( 15f)  is  ob¬ 
tained  by  matching  the  solution  in  the  Hartmann  layer 
between  the  free  shear  layer  and  the  insulating  wall. '  ^  Physi¬ 
cally,  this  Hartmann  layer  must  match  the  free 

shear  layer  velocity  Uy,  o  at  r  =  0,  and  must  satisfy  the  no¬ 
slip  condition  at  the  insulating  wall.  This  jump  in  velocity 
involves  a  very  large  gradient  in  the  viscous  shear  stress  in¬ 
side  this  Hartmann  layer,  i.e.,  is  large.  The  term  t„ 

represents  the  x  component  of  the  surface  force  per  unit  area 
along  the  x-y  plane.  An  0(M ' '^)  tangential  current  density 
jy  inside  this  Hartmann  layer  provides  the  electromagnetic 
body  force  to  balance  this  large  shear  stress  gradient.  The 
total  Hartmann  layer  current  in  the;’  direction  at  each  point 
along  the  insulating  wail  is  proportional  to  the  free  shear 
layer  velocity  «y,,o  evaluated  at  f  =  0.  Since  this  velocity 
varies  along  the  insulator,  the  total  current  inside  the  Hart¬ 
mann  layer  must  vary  as  well.  An  increase  in  current  must  be 
drawn  from  the  free  shear  layer,  while  a  decrease  implies  a 
current  into  the  free  shear  layer  from  the  Hartmann  layer.  In 
the  boundary  condition  ( 1 5f ) ,  the  term  d  ^  ^/d^  ^  represents 
the  tangential  variation  of  «/,.o  at  t  =  0  since  «yi.o  =  M>/ 
d^.  The  term  d^/dt  represents  the  electric  current  from  the 
Hartmann  layer  into  the  free  shear  layer  or  vice  versa,  since 

y,  =  - 

By  using  the  techniques  of  complex  Fourier  transform, 
the  convolution  theorem  can  be  used  to  obtain  ‘t>(t,f): 


t  =  j]  +  a  sin  0, 
and 


(I4b) 


sin^>0  for /I,  sign(sin  0)  =  -|- 1.  (14c) 


Here  ‘I>(/,f)  is  treated  as  a  perturbing  term  which  ap¬ 
proaches  zero  as  -  00  [see  boundary  value  Eq.  (9a)  ]. 
Also  approaches  zero  at  the  moving  perfect  conduc¬ 
tor.  The  partial  differential  Eq.  (13b)  in  terms  of  <t>(r,f) 
with  its  five  boundary  conditions  in  terms  of  t  and  ^  is 


d-<t> 

dt~ 


d*d> 

dC* 


for  0</<sec  d. 


—  00  <C<  <»• 

(15a) 


‘l>=r  " 

J  ao  IT 

=  r  'V(^*)^('.f-?*K*.  (16a) 

Jo  —  IT 


where 


>  d^<l>  d(l>  l< 

/(f)=-^(0.f)+^(0,f)  =  | 


0  for  —  00  <0, 

/(f)  forO<f<oo. 

(16b) 


In  the  previous  work'^  the  Green  function  was  evaluat¬ 
ed  to  have  the  form 
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—  erf[Jf(2  sec  0  —  /)“  '^"1} 

+  i;r-'^Hr'^=exp(-if^r-') 


—  (2  sec  0  —  t)'^-  exp{  —  ’(2  sec  0  —  t)  ']}. 

(17) 

This  solution  satisfies  the  governing  equation  and  all  bound¬ 
ary  conditions  except  the  condition  ( 15f). 

Remembering  that  can  be  expressed  by  the  con¬ 

volution  integral  [Eq.  (16a)]  and  must  also  satisfy  the 
boundary  condition  at  the  fixed  perfect  conductor 


‘b  =  —  (<^„  —  a  cos  0)  at  t  =  0,  for  0 < ^ <  oo. 


(18) 


the  following  relationship  must  hold: 

-(.S^,-acos0) 


L 


forO<f<oo.  (19) 

In  this  equation,  we  change  coordinates  to  reduce  the  inte¬ 
gral  equation  to  a  form  that  is  independent  of  sec  0  and  of 
(<^i,  -  a  cos  0). 


i 


/’(Z*){|Z-Z*|erfc(|Z-Z*|) 

-7r-''=exp[  -  (Z-Z*)=]}dZ*  =  1 
for  0  <  Z  <  00 , 
where 

;  =  2(2 sec 0)''- Z,  ;*  =  2(2 sec 0)''-Z*, 
rff*  =  2(2sec0)'^-dZ*,  0<f,  Z<oo, 


(20a) 

(20b) 

(20c) 


/(^)=/[2(2sec^)"=Zj= 

2  sec  0 

(20d) 

This  integral  equation  is  solved  numerically  in  detail  here  in 
order  to  determine  the  velocity  profile  directly  in  the  free 
shear  layer.  We  did  not  solve  it  in  the  previous  paper.'* 

The  velocity  of  the  free  shear  layer  is  represented  to  zero 
order  by  [see  Eq.  ( 12b)  ] 

u  =  M'‘-  (21a) 

where  the  stretched  coordinate  f  is  defined  as 

;  =  M''\^  +  aco%0).  (21b) 

Here  Uj  |,„  is  related  to  the  zero-order  potential  function  by 
0<P 


(U) 


forO</<sec0,  —  oo<^<oo. 
The  function  „  can  be  represented  as 


(22) 


«y>.o=  ['  (23a) 

Jtt  dg 


where 
dG 


(/.b)  =  i{erf(].C/  "=) 

-erf[l.C(2sec6>-/)  ''']} 


(23b) 


was  obtained  from  G(t,g)  by  partial  differentiation  [see  Eq. 
( 17)  ].  By  using  the  identities 

erf(ac)  =  sgn(jc)erf(  |x| ),  erf(x)  =  1  -  erfc(x). 


we  obtain  the  expression 
4  Jo 

x{erfcij|f  —  ^*|(2sec0  —  r)-’'’^] 
-erfc((|f-^*|t-''^)K*,  (24a) 

where 

f(.g)=  —  U<J„  — flcos0)cos0F[lf(2sec^)“''-].  (24b) 
Substituting 

Z={g{2  sec  0)  -'^-*-*g  =  2(2  sec  0)''-  Z,  (25a) 

Z  *  =  It"  *(2  sec  0)  -''-<-*g  *  =  2(2  sec  0)''~Z  *,  (25b) 

T=\  cos  0t->^t  =  2  sec  0T.  (25c) 

into  Eq.  (24a)  results  in  the  expression  for  U/i.o: 

“/i.o  =  (<!>{)  — a  cos  0) (cos  0)''-  U{ T,Z),  (26a) 

where 

U=(2)-^>'-r  F(Z*)sgn{Z-Z*) 

Jo 

x{erfc(|Z-Z*|7’-'^^)  -erfc 
X[|Z-Z*|(1 -D-'/^lVZ*.  (26b) 


It  remains  to  solve  the  integral  equation  for  F{Z*)  for 
0<Z  *  <  00  and  to  substitute  this  solution  into  the  integral 
expression  for  U{T,Z).  The  universal  velocity  function 
U{  T,Z)  for  0<  r<0.5  and  —  oo  < Z<  oo  is  independent  of 
all  parameters  in  the  problem.  For  a  particular  set  of  values 
for  0  and  ^o,  expression  (26a)  scales  this  universal  velocity 
function  into  the  0(.M  '^^)  free  shear  layer  velocity  Uy,  q. 

The  solution  of  the  integral  equation  involves  a  Dirac 
delta  function  at  Z  *  =  0.  The  Hartmann  layers  on  the  insu¬ 
lator  and  perfect  conductor  intersect  to  form  an 
corner  region  aty  =  0,  z=  —a,  as 
shown  in  Fig.  4.  As  noted  in  part  a  of  this  figure,  the  potential 
changes  from  0  to  ((i„  —  o  cos  0)  over  a  very  small  distance 
Ay  =  atz=  —o.  This  is  the  same  as  the  jump  in 

<f>  across  the  free  shear  layer  at  t  =  0,  but  in  the  free  shear 
layer  this  potential  difference  is  balanced  by  the  large  in¬ 
duced  electric  field  (wy,  gx)  Xb.  Inside  the  Hartmann  layer 
adjacent  to  the  insulator,  u  must  decrease  from  «y,,o  to  0,  so 
that  locally  the  large  potential  gradient  is  not  balanced  by 
the  induced  electric  field.  This  leads  to  a  large  electric  cur¬ 
rent  leaving  the  fixed  perfect  conductor  through  the 
0(  A/  “ ' )  X  0(  A/  “ ' )  corner  region.  Part  of  this  current  con¬ 
tinues  in  the  y  direct  to  become  the  jy  in  the  insulating  wall 
Hartmann  layer  which  balances  the  large  viscous  shear  gra¬ 
dient  here.  The  rest  of  this  current  enters  the  free  shear  layer 
directly.  Since  the  corner  region  has  0{M  ~ ')  dimensions, 
while  the  free  shear  layer  has  a  much  larger  thick¬ 

ness,  the  free  shear  layer  sees  the  current  from  the  corner 
region  as  a  point  source  of  electric  current  at  r  =  0  and  f  =  0. 
Since  Fis  the  rescaled  value  of 

^  (27) 

dg-  dt 
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0(M-'l  X  0(M-'l 
comer  region 


(b)  (c) 


FIG.  4.  Hartmann  layers  whrch  separate  the  free  shear  layer /I  from  the 
walls  at  the  corner,  the  0(  A/  '}xO(M  ')  corner  region  and  the  electric 
current  flowing  through  the  corner  region  into  the  free  shear  layer. 


at  t  =  0,  and  y,,  =  —  d<i>/dt  is  the  current  which  sees  a 
source  at  Z  =  0,  then 

^{Z*)  =/„<5(Z*-f) +F(Z*),  (28) 

where  0<e<l  is  included  to  avoid  any  ambiguity  in  the 
integral  of  S  from  Z  ♦  =  0  to  oo .  Here  /„  represents  the  mag¬ 
nitude  of  the  electric  current  flowing  from  the  corner  region 
to  the  free  shear  layer  and  is  determined  by  the  condition 
that  F(  Z  * )  is  bounded  at  Z  *  =  0.  In  other  words,  /„  '*  cho¬ 
sen  to  remove  all  singular  behavior  from  the  modified  inte¬ 
gral  equation 


*){|Z-Z*|erfc(|Z- 


Z*\) 


-7r-'^’exp[  -  (Z-Z*)-]}ciZ* 


=  l-t-/„[;r  '^*exp{  —  Z‘)  —  Zerfc(Z)].  (29) 


This  integral  equation  is  solved  numerically  using  the 
following  steps.  The  integral  from  Z  ♦  =  0  to  ZN  is  approxi¬ 
mated  by  a  trapezoidal  rule  and  with  a  step  size  AZ  *  =  H. 
For  the  integral  from  Z*  =  ZN  to  oo,  we  assume  that 
F(Z*)  =  F(ZN)  and  carry  out  the  integration  explicitly. 
The  integral  from  ZN  to  oo  is  then  added  to  the  coefficient  of 
F{ZN)  from  the  trapezoidal  rule  quadrature  formula.  The 
integral  equation  is  evaluated  at  Z  =  iH,  for  /  =  0  to 
n  =  ZN/H,  i.e.,  at  Z  =  0,  H,  2H,...ZN.  This  gives  n  -f  I  si- 
multajneous  linear  algebraic  equations  for  the  unknowns 
Fj  =  F(jH)  for j  =  0ton.  This  set  of  simultaneous  equations 


FIG.  5.  Local  (r,^)  coordinates  for  each  layer/1  and /2.  ( 1 )  Origin  is  at  the 
comer  of  conductor  and  insulator,  ( 2 )  r  is  parallel  to  the  magnetic  field;  ( 3 ) 
at  corner,  perfect  conductor  lies  m  ^  >  0  and  insulator  lies  in  ^  <  0,  and  (4)  ^ 
is  the  regular,  unstretched  coordinate. 

is  solved  using  Gauss  eliminadon  with  partial  pivoting.  The 
value  of  /„  is  adjusted  until  Fy  =  F(0)  is  clearly  the  finite 
continuation  of  F,  =  F{H)  and  Fj  =  F(2H).  The  values  of 
Fy  for  two  successive  values  of  Iq  were  used  to  predict  the 
next  guess  for  /„.  This  process  converges  quickly  to 
/q  =  —  1 .  This  test  is  quite  sensitive  because  Z,,  =  —  0.98 
gives  a  very  large  positive  value  of Fy  and  /y  =  —  1 .02  gives  a 
very  large  negative  value  of  Fy.  For  /„=  —1.00, 
F„=-1.48,  F,  =  -1.51,  and  F2=-1.54,  with 
H  =  0.03.  The  values  of  H,ZN,  and  n  were  varied  over  a  wide 
range  of  values  to  test  the  accuracy  of  the  trapezoidal  rule 
quadrature.  The  results  for  H  =  0.03,  ZN  =  1 .2,  and  «  =  40 
were  the  results  used  to  compute  U,  but  were  essentially 
identical  to  the  values  for  the  other  similar  sets  of  H,  ZN,  and 
n. 

The  expression  for  C/(  T,Z)  is  evaluated  in  a  comps  table 
series  of  steps.  We  substitute  F{Z*)=  —  S(Z*  —  e) 
+  F(Z*).  We  use  a  trapezoidal  rule  quadrature  for 
0<Z *<ZA^.  WetakeF(Z *)  =  F„  forZ *  >  ZA'and integrate 
for  this  range  explicitly. 

IV.  GENERALIZED  FREE  SHEAR  LAYER  SOLUTIONS 
A.  Free  shear  layer  f1 

We  have  defined  to  equal  —  a  cos  0)  in  free  shear 
layer  fl  (see  Figs.  5  and  6).  The  change  in  potential  across 
the  layer  is 


FIG.  6  Free  shear  layer /I. 
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^(f>(,t)  =  tx))  —  —  m),  (30a) 

where  the  stretched  coordinate  Z  is  defined  as 

^  =  2A/-''=(2sec(9)''-Z,  (30b) 

2  =  \M''-{litce)-'%  (30c) 

where  ( t,^)  are  local  Cartesian  coordinates  with  the  origin  at 
the  comer.  Therefore, 

A^(/)  =^,(/cos0- 1),  (31) 

so  that  A^(r)  varies  from  —  at  r  =  0  to  0  at  /  =  sec  'fhe 
velocity  in  the  free  shear  layer  /I  is  expressed  as 

u  =  3/''V.(cos0)'/=f/(r,Z),  (32) 

where  t  has  been  defined  as  t  =  T{1  sec  6)  in  Eq.  (25c). 
r  =  0  at  the  corner  of  the  insulator  and  perfect  conductor 
and  T  =  0.5  at  the  other  conductor.  The  range  of  T  is 
0<7’<0.5,  which  is  inconvenient,  but  it  eliminates  several 
square  roots  of  two  from  the  equations. 

When  Z—  00  the  stretched  coordinate  approaches  the 
core  between  two  perfect  conductors.  When  Z—  —  oo,  the 
coordinate  approaches  the  core  between  the  conductor  and 
the  insulator.  U(T,2)  is  independent  of  0,  and  a. 

B.  Free  shear  layer  f2 

In  free  shear  layer /2,  (f>2  was  defined  to  equal 
<f>2  =  (^<.  +  a  cos  9  —  sin  6) 

(see  Figs.  5-7).  The  jump  in  potential  across  the  layer  is 
A^(/)  =  ^(Z-»  00 )  —  ^(Z-*  —  00 ),  (33a) 

where  the  stretched  coordinate  Z  is  defined  as 

|■  =  2A/-''^(2sec(9)''^Z,  (33b) 

2='2M''H2sec0)-'%  (33c) 

so  that 

A(^(r)  =  -<A2(tcos0- 1).  (33d) 

The  velocity  in  the  free  shear  layer /2  is  expressed  as 

u=  -M''%icosd)'>W(T,2).  (34) 

The  definitions  for  T,2,  and  U  are  the  same  as  those  given 
previously. 

C.  Generalized  expression  for  velocity  in  free  shear 
layers 

The  change  in  potential  A^(r)  across  the  generalized 
layer  is  (see  Fig.  8) 

A^(r)  =  ^(Z-»  00 )  —  ^(Z->  —  00 ),  (35) 
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and  hMulator 
at  any  angle 


J  >  0;  parfact  conductor 
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FIG.  8.  Generalized  notation  for  free  shear  layer. 


where  Z  is  the  stretched  coordinate.  Then 

A(^U)  =<^i,  +  i<f>,  -<l>i 

=  ('!>, (36a) 
and  the  velocity  u  in  the  free  shear  layer  is 

u  =  ''^U{ T,2),  (36b) 

where 

t=^2t,T*-*T=t/2t„  (36c) 

and  the  stretched  coordinate  Z  is  defined  as 
i  =  2M-  )'''2^2  =  J  M  '^^(2f, ) " 

(36d) 

As  Z-*oo,  the  stretched  coordinate  approaches  the  core 
between  the  two  perfect  conductors,  and  as  Z-*  —  oo  it  ap¬ 
proaches  the  core  between  the  conductor  and  the  insulator. 

The  function  U[T,2)  was  defined  as  before.  The  results 
for  U are  presented  in  Fig.  9.  Since  inside  the  free  shear  layer 

u  =  ^,  ^  =  (37) 

and  since  ^  =  0  at  /  =  0,  for  f  >0,  then  k  =  0  at  f  =  0,  for 
f  >0.  As  7’-«0,  the  velocity  profile  becomes  skewed  to  nega¬ 
tive  Z  with  K— 0  for  Z>0.  The  volume  flux  inside  the  free 
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FIG.  7.  Free  sheer  layer  /2. 


FIG.  9.  Rescaled  free  shear  layer  velocity  profiles  (/( T,Z). 
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shear  layer  is  proportional  to  (sec  d  —  t),  so  it  is  0  at 
t  =  sec  ^  or  T=  0.5,  and  maximum  at  t=T=0  where 
u  =  0  for  f  >0.  As  T  increases  from  0,  the  profile  settles  to¬ 
ward  a  nearly  symmetric  one  with  its  center  at  Z  =  —  0.3. 

V.  VISCOUS  DISSIPATION  AND  JOULEAN  POWER 
LOSSES 

A.  Free  shear  layers 

The  nondimensional  viscous  dissipation  in  the  total 
channel  can  be  found  by  solving  the  following  integral  in  the 
different  regions  (see  Fig.  10): 


(38) 

The  contributions  of  the  core  regions  are  of  order  0(1),  the 
contributions  of  the  Hartmann  layers  between  the  cores  and 
the  walls  are  of  0(  3/) ,  the  contributions  of  the  two  free  shear 
layers  are  of  order  0(  M ,  ^nd  the  contributions  of  the  two 
Hartmann  layers  between  the  free  shear  layers  and  the  insu¬ 
lating  walls  are  of 

The  nondimensional  viscous  dissipation  in  free  shear 
layer  /I  is  expressed  as 

/■w,  (39a, 

where 

«^i,o(r,f)  =  (4 -a cos 0) (cos ^) '^=6/(7', Z)  (39b) 
and 

=  (2)  -  ’^'(4  -  a  cos  0)cos  6  T,Z). 

dg  aZ 

(39c) 

Since  df  =  2(2  sec  0) dZ  and  dt-=2  sec  6  dT,  the  vis¬ 
cous  dissipation  in  /I  can  be  e^'pressed  as 


FIG  10.  Subregions  m  magnetohydrodynamic  channel  flow  at  high  exter¬ 
nal  magnetic  fields 
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Ptm/i  =  -a  COS  0f(cos  (40a) 

where  i 

We  introduce  Eq.  (26b)  for  U(T,Z),  and  expression  (28) 
for  F{Z*).  We  use  the  identities 

A(i2|)=sgn(Z),  4:tsgn(Z)]  =25(Z).  (41) 

aZ  aZ 

Finally  we  introduce  the  numerical  results  for  F(Z*)  and 
evaluate  the  integral  with  respect  to  Z  *  with  the  same  trape¬ 
zoidal  rule  quadrature  to  obtain  a  value  for  P^^j-  is  a 
universal  constant  for  all  free  shear  layers  of  this  type  and 
equals  0.09 1 8.  The  dimensionless  viscous  dissipation  for  free 
shear  layer  f\  is  written  as 

Pt,,j  I  =  (cos  0)  "-(P„/),  (42) 

where  i^,  =  —  a  cos  0). 

The  dimensionless  viscous  dissipation  for  free  shear  lay¬ 
er  fl  is  derived  from  analogous  arguments  as 

Pr,.fi  =  M  (cos  0)''-{P,f),  (43) 

where  =  (i^,,  +  a  cos  0  —  sin  0). 

Also  the  viscous  dissipation  for  a  generalized  free  shear 
layer  using  the  generalized  formulation  for  the  free  shear 
layers  can  be  expressed  as  [Eqs.  (36)  ] 

Pj,.f  =  M^'H4>,-<l>S'U,)-''\P^r)>  (44) 

where  P^^f  is  again  the  universal  constant  0.0918. 

B.  Viscous  dissipation  in  the  Hartmann  layers 

Let  P|■^,H  1  be  the  dimensionless  viscous  dissipation  for 
the  Hartmann  layer  between  the  free  shear  layer  /I  and  the 
insulator  at  z  =  —  a  and  let  be  the  dimensionless  vis¬ 
cous  dissipation  for  the  Hartmann  layer  between  the  free 


FIG.  11.  Hartmanniayeradjacent  to  free  shear  layer/!.  The  0(  A/ free 
shear  layer  velocity  u^,  „  is  zero  at  r  =  sec  0  for  all  ^  and  at  /  =  0  for  0. 
There  are  no  0(  A/ ' ' ' )  velocities  in  the  Hartmann  layers  between  the  free 
shear  layer  and  the  perfect  conductors,  so  that  =  0(A/''-)intheselay- 

ers. 
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shear  layer /2  and  the  insulator  at  z  =  a. 

In  a  previous  paper,'*  we  found  the  zero-order  velocity 
in  the  Hartmann  layer  /f  (see  Fig.  1 1 )  to  have  the  form 

u  =  M  (0,  —  Y  sin  d) 

X[1 -exp(  -  |sin0|=)],  (45) 

where 

Y=M''^y  and  E  =  M{z  +  a)  (46) 

are  the  stretched  tangential  and  normal  coordinates  for  this 
Hartmann  layer,  respectively.  Here  f  =  —  F  sin  0  in  the 
free  shear  layer  at  /  =  0.  In  the  Hartmann  layer 

-^  =  —  Fsin  0)|sin0| 

dz 

Xexp(  -  |sin  01=),  (47) 

while  0u/(?y  is  only  0{M).  We  introduce  the  result  Eq.  (47) 
into  expression  (38)  for  the  dimensionless  viscous  dissipa¬ 
tion,  where 

dy  =  M-''-dY,  dz  =  M-'dE,  (48) 

for  this  Hartmann  layer.  The  ranges  of  integration  for  both  Y 
and  H  are  from  0  to  oo .  The  integral  with  respect  to  H  can  be 
carried  out  explicitly  since  it  only  involves 
exp(  —  2|sin  6\B).  The  result  is  an  expression  for  the  di¬ 
mensionless  viscous  dissipation  in  this  Hartmann  layer, 

P,„„^  =(A/''-isin0|  r  rsin0)]-rfF. 

Jo 

(49) 

Substituting  -  Fsin^,  or  £?F=  —csc6d^,  noting 
that  sin  0  is  po.sitive  for  this  Hartmann  layer,  and  reversing 
the  limits  of  integration  gives 

(50) 

Note  that  the  expression  is  independent  of  the  angle  of  the 
insulator  relative  to  the  free  shear  layer.  The  integral  in¬ 
volves  only  the  free  shear  layer  solution  which  depends  on 
only  two  parameters:  the  jump  in  0,  namely 
=  a  cos  0),  and  the  length  of  the  layer  t,  =  sec  0. 
Changing  the  variables  from  u,  i.o  to  (/ and  from  y  to  Z  gives 


P„in  =M''-<l,]{cos6)'^HP,,„), 

(51a) 

where 

=  (2)''-  f’  (C/(0,Z)]-r/Z. 

J  ~  Ob 

(51b) 

Pi,„  was  calculated  numerically  from  Eq.  (26b)  and  the  val¬ 
ue  is  0.197  26,  which  is  more  than  double  the  value  for  the 
free  shear  layer. 

(7(0, Z)  has  a  discontinuity  or  step  at  Z  =  0  from  rough¬ 
ly  —  0.7  to  0,  as  shown  in  Fig.  12.  The  corner  region  which  is 
0(M  ')X0(M  ')  must  match  this  discontinuity.  Hunt 
and  Stewartson  give  the  solution  inside  the  corner  region. 
In  this  region  u  —  0{M'''),du/dy  =  and  du/ 

dz=  so  that 

However,  dy  =  0{M  ' )  and  dz  =  0(M  ' )  in  expression 
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FIG  12.  Sketch  of  t/(r,Z)  at  7'=  0. 


(38)  for  the  viscous  dissipation,  so  that  the  contribution  to 

from  the  corner  region  is  only  O(M').  This  contribution 
represents  a  small  correction  to  the  viscous  dissipation  in  the 
Hartmann  layer  between  the  insulator  and  free  shear  layer 

/I. 

The  only  0(M  contributions  to  the  viscous  dissipa¬ 
tion  arise  from  the  two  free  shear  layers  and  the 
0{M~')xO{M~''~)  Hartmann  layers  separating  the  free 
shear  layers  from  the  adjacent  insulators.  Summing  both  free 
shear  layers  and  their  Hartmann  layers  gives  the  following 
expression  for  the  viscous  dissipation 

■\-6\){cosdy'HP,ft-P„H)>  (53) 

where  P„;  =  0.0918  and  P„„  =0.197  26. 

The  coefficient  (F^y  -b  P^h  )  =  0.2891  of  the 
dimensionless  viscous  dissipation  is  relatively  small.  With 
the  assumption  that  1,  we  have  so  far  assumed  that  the 
0{M)  contributions  to  the  viscous  dissipation  are  negligible 
compared  to  the  leading  term  in  the  asymptotic 

expansion.  In  actual  homopolar  devices,  A/  has  a  moderately 
large  value,  for  example,  50,  so  that  M is  not  much  larger 
than  M.  If  the  coefficient  of  the  0{M)  term  is  larger  than 
(F^,y  -I-  Pi,n ),  then  the  0(M)  term  might  be  as  large  as  the 
0{M  ’'-)  term  for  actual  values  of  M,  such  as  50.  There  are 
four  contributions  to  the  0(M)  dimensionless  viscous  dissi¬ 
pation: 

(1)  The  0(M)  correction  or  perturbation  in  the  free 
shear  layers. 

(2)  The  0{M)  correction  or  perturbation  in  the  Hart¬ 
mann  layers  separating  the  free  shear  layers  from  the  insula¬ 
tors. 

( 3 )  The  viscous  dissipation  in  the  0(  A/  ~  ‘ )  X  0(  A/  “ ' ) 
corner  regions  atji  =  0,  z  =  ±a,  associated  with  the  jump  in 
(/(0,Z)  from  —  0.7  to  0. 

(4)  The  viscous  dissipation  in  the  Hartmann  layers 
between  the  core  regions  and  the  walls  where  there  is  a  jump 
in  the  0(  1 )  velocity  across  the  Hartmann  layer. 

The  first  three  contributions  to  the  0{M)  dissipation 
are  all  proportional  to  the  free  shear  layer  solution.  We  will 
assume  that  this  proportionality  leads  to  coefficients  which 
are  comparable  to  (F^,y  -I-  ),  so  that  these  contributions 

are  indeed  0(  A/  *  ' '  - )  smaller.  On  the  other  hand,  the  Hart- 
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mann  layers  adjacent  to  the  core  regions  involve  a  contribu¬ 
tion  that  is  related  to  the  potentials  in  the  perfect  conductors 
and  not  to  those  in  the  free  shear  layers.  In  particular,  the 
viscous  dissipations  in  the  Hartmann  layers  adjacent  to  the 
two  perfect  conductors  increase  as  the  aspect  ratio  of  the 
channel,  2a,  increases.  Real  current  collectors  have  large 
values  of  2a,  for  example,  25,  so  that  the  <?(A/)  Hartmann 
layer  contribution  may  be  larger  than  the  free  shear 

layer  contribution. 

The  Hartmann  layers  adjacent  to  the  cores  are  num¬ 
bered  A  1  to  A  6  in  Fig.  10.  The  core  solutions  were  presented 
in  a  previous  paper. Hartmann  layers  A  1,  /i  3,  and  h  4  in¬ 
volve  no  jumps  in  the  0{  1 )  velocity  u,  and  their  viscous 
dissipation  is  at  most  0(  1 ) .  For  Hartmann  layers  A  2,  A  5,  and 
A  6,  the  velocity  varies  exponentially  from  the  core  value  to 
the  wall  value  ( 0  or  1 ) .  The  velocities  inside  these  Hartmann 
layers  are 

«A6  =  1 -exp[ -4/|sin0|(z  +  a)],  (54a) 

“*2  =Wc2(>'=  1){1  -exp[il/|cos5|(y-  1)]}+  1, 

(5'-  ;) 

“/.5  =  «c2(.V  =  0){l  -exp(  -A;|cos0|>>]}.  (54c) 

In  core  c2,  (j)  is  determined  by  the  potentials  of  the  perfect 
conductors  and  u  =  dij)/d^.  Integrating  the  viscous  dissipa¬ 
tion  per  unit  volume  over  each  of  the  Hartmann  layers  A  6, 


A  2,  and  A  5  gives 

P^h<,=\Msme,  (55a) 

(55b) 

where  again  <f)f  =  <f>o  -  a  cos  6,  and 

^2  =  ^0  +  a  cos  0- sin  0  =  -1- A^,  (55c) 

A^  =  2fl  cos  0  — sin  0>O,  (55d) 

for  0<tan  6<,2a.  Therefore,  the  viscous  dissipation  in  the 
free  shear  layers  and  adjacent  Hartmann  layers  /f  is  propor¬ 
tional  to 

2<l>]  +  2<i,(A^)  -I-  (A^)',  (56) 


while  the  viscous  dissipation  in  the  Hartmann  layers  adja¬ 
cent  to  the  core  cl  is  proportional  to 

(l>\-<f>\  =  +  3^,(A^)  +  (A^)'],  (57) 

The  quantity  (A^)  can  be  moderately  large.  For  example, 
for  2a  =  25  and  9  =  45*,  ( A<^)  =  17,  which  can  make  the 
Hartmann  layer  dissipation  comparable  to  that  of  the  free 
shear  layers  for  modest  values  of  M.  Therefore,  the  best  ap¬ 
proximation  to  the  viscous  dissipation  would  be  the  sum  of 
the  contribution  from  the  free  shear  layers  and 

Hartmann  layers  H  and  of  the  0(M)  contribution  from  the 
Hartmann  layers  A  6,  A  2,  and  A  5. 

C.  Joulean  power  lOMes 

In  addition  to  viscous  dissipation,  there  are  Joulean 
power  losses  due  to  the  flow  of  electric  current  through  the 
electrically  resistive  liquid  metal.  The  Joulean  power  losses 
are  given  by 

Pv  =  ^  =  M^[  f  (//  ^n)dy  dz.  (58) 
II/Uq  J-aJo 
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In  the  core  regions  cl  and  c3,y,  =  0(M~');  in  the  core 
region  c2,  y,  =0(1);  and  in  all  three  core  regions, 
yj  =  0(.M~^).  Therefore,  the  contributions  to  the  Joulean 
power  losses  are  O(M^)  for  the  core  region  c2  and  0(  1 )  for 
the  core  regions  cl  and  c3.  In  a  previous  paper,  we  present¬ 
ed  the  expression  for  the  Joulean  losses  in  the  core  region  c2, 
namely 

Pr^2  =  i  V "(cos 9) (59) 

Equation  (59)  assumes  that  y,,  =  —  cos0(^o  +  l’)  every¬ 
where  in  the  parallelogram  defined  by  (tan  0)|'<i/< 
(tan  9)^  -I-  sec  9,  for  —  a  cos  9<i<a  cos  9  —  sin  9,  and 
that  j  =  0  everywhere  in  the  two  triangles  for  !■<  —  a  cos  9 
and  for  cos  9  —  sin  9.  In  reality,  the  0(  1 )  current  den¬ 
sity  deviates  from  these  values  inside  the  free  shear  layers  /I 
and /2  and  in  the  Hartmann  layers  A  2  and  A  5,  so  that  these 
viscous  layers  contribute  corrections,  either  additions  or 
subtractions,  to  the  Joulean  power  losses  based  on  the  core 
solutions. 

The  current  lines  do  not  actually  end  exactly  at 
^  =  —  a  cos  0  and  &X^  =  acos9  —  sin  9.  Instead,  the  cur¬ 
rent  lines  fringe  an  0{M  distance  beyond  these  lines 
inside  the  free  shear  layers.  Thus  the  cross-sectional  area  of 
the  effective  conductor  between  the  two  electrodes  is  the 
0{  1 )  area  of  the  core  region  c2  plus  the  0(M  "  '^^)  areas  of 
the  two  free  shear  layers  /I  and  /2.  The  addi¬ 

tional  current  inside  the  two  free  shear  layers  contributes  an 
additional  Joulean  power  loss.  The  steps  to  deter¬ 
mine  the  free  shear  layer  contribution  exactly  parallel  the 
steps  to  determine  the  viscous  dissipation  due  to 

the  free  shear  layers.  The  total  addition  due  to  both  free 
shear  layers  is 

Prif  =  0.m2M^'^ (cos  9)' (60) 

The  additional  Joulean  power  losses  due  to  the  extra  fringing 
current  inside  the  free  shear  layers  is  more  than  twice  the 
viscous  dissipation  in  the  free  shear  layers  and  associated 
Hartmann  layers,  as  given  by  Eq.  (53). 

In  the  Hartmann  layers  A  2  and  A  5,  the  current  lines 
from  the  core  region  c2  must  bend  in  order  to  enter  the  per¬ 
fect  conductors  at  right  angles.  This  bend  in  the  current  lines 
makes  the  current  density  in  the  Hartmann  layer  less  than 
that  in  the  core.  Since  Eq.  (59)  assumes  that  the  core  solu¬ 
tion  applies  all  the  way  to  the  perfect  conductors,  it  slightly 
overestimates  the  Joulean  power  losses  very  near  the  con¬ 
ductors.  The  0(M)  correction  due  to  both  the  Hartmann 
layers  A  2  and  A  5  is  given  by 

Fjj*  =  —  Af(sin  9)^  cos  9(4)1  —  <l>\  )•  (61 ) 

The  best  estimate  of  the  Joulean  power  losses  is  given  by  the 
sumofEqs.  (59)-(61). 

VI.  CONCLUSIONS 

In  a  rectangular  duct  having  insulating  sides,  a  moving 
perfectly  conducting  top  wall,  a  stationary  perfectly  con¬ 
ducting  bottom  wall,  and  a  strong,  skewed  magnetic  field, 
there  are  free  shear  or  interior  layers  lying  along  the  two 
magnetic  field  lines  through  opposite  corners  of  the  cross 
section.  These  free  shear  layers  separate  pairs  of  inviscid  core 
regions  and  must  match  different  electrical  potentials  in  the 
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adjacent  cores.  These  layers  involve  very  large,  0(M  di¬ 
mensionless  velocities,  while  the  core  velocities  are  only 
0(  1 ).  These  high  velocity  sheet  jets  can  be  in  either  the  plus 
or  the  minus  x  direction,  i.e.,  with  or  opposite  to  the  moving 
wall’s  direction.  The  direction  of  flow  in  each  free  shear  layer 
depends  on  the  sign  of  the  electrical  potential  at  the  point 
where  the  layer  intersects  the  moving  conductor,  that  is,  on 
the  signs  of  and  <t>2-  These  values  depend  on  the  direction 
of  the  load  current:  if  there  is  a  significant  load  current  from 
the  stator  to  the  rotor,  then  both  (f>^  and  ^2  are  negative;  if 
there  is  a  significant  current  from  the  rotor  to  the  stator,  then 
both  ij)  I  and  (4,  are  positive;  if  there  is  no  load  current,  then 
is  negative  and  ^2  is  positive. 

The  boundary  value  problem  for  the  free  shear  layers 
was  reduced  to  an  integral  equation  in  a  previous  paper.*’ 
The  numerical  solution  of  this  integral  equation  is  presented 
here.  The  velocity  profiles  in  the  free  shear  layers  are  also 
presented  here.  These  profiles  are  reduced  to  a  universal 
function  U{T,Z)  for  all  similar  shear  layers.  The  scaling  for 
any  particular  shear  layer  depends  on  only  two  parameters: 
the  jump  in  the  electrical  potential  across  the  layer  at  the 
corner  where  the  insulator  and  conductor  meet,  i.e.,  and 
^2  for  free  shear  layers  /I  and  fl  here,  and  the  length  of  the 
free  shear  layer  along  the  magnetic  field  line,  i.e.,  sec  0  here. 

Since  there  are  large  velocities  inside  very  thin  regions, 
the  free  shear  layers  and  adjacent  Hartmann  layers  involve 
large  viscous  dissipation.  Their  contributions  to  the  dimen¬ 
sionless  viscous  dissipation  are  0{M  '^’j,  while  those  of  the 
cores  are  only  0(  1 ).  We  have  also  computed  the  0(A/)  vis¬ 
cous  dissipation  in  the  Hartmann  layers  between  the  cores 
and  the  walls.  Finally  we  determined  the  contributions  to  the 
Joulean  power  losses  from  the  core  region  cl,  from  the  free 
shear  layers/ 1  and  fl,  and  from  the  Hartmann  layers  h  1  and 
h5. 

Here  we  assumed  that  the  flow  is  laminar.  Very  strong 
magnetic  fields  would  completely  eliminate  turbulence  and 
instability  in  the  laminar  flow.  However,  in  actual  homopo- 
lar  devices,  the  magnetic  field  is  strong  enough  to  make  the 
turbulence  two-dimensional,  but  is  not  strong  enough  to 
eliminate  turbulence  completely  or  to  eliminate  instabilities 
which  lead  to  vortices  aligned  with  the  magnetic  field.  The 
two-dimensional  turbulence  or  aligned  vortices  would  be 
stnmgest  m  the  free  shear  layer  because  of  the  large  velocity 
gradients  here.  The  effects  of  this  turbulence  or  voriicity 
would  be  to  increase  the  thickness  of  the  free  shear  layer,  to 
reduce  the  maximum  free  shear  layer  velocity  proportion¬ 
ately,  and  to  reduce  the  viscous  dissipation  considerably. 
Two-dimensional  turbulence  or  aligned  vortices  in  a  strong 
magnetic  field  involve  relatively  little  dissipation,  unlike 
three-dimensional  turbulence.  Therefore,  the  present  ex¬ 
pressions  for  the  viscous  dissipation  should  provide  conser- 
vati'.e  overestimations. 
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APPENDIX:  HUNT  AND  STEWARTSON’S  INTEGRAL 
EQUATION  FOR  A  SIMILAR  FREE  SHEAR  LAYER 


The  integral  equation  of  interest  in  the  present  paper  is  a 
Fredholm  integral  equation  of  the  first  kind 

rF(Z*){|Z-Z*|erfc(|Z-Z*|) 

Jo 

-5r-''-exp[  -  (Z-Z*)-]VZ*=  1.  (Al) 

Hunt  and  Stewartson’s'*’  integral  equation  has  the  form 


g{p)  =  (In) 


g(0exp[ 

0 


k(t-p)-]dt.  (Al) 


Equation  (Al )  can  be  reduced  to  Eq.  (A2)  by  three  steps: 
integration  by  parts,  differentiation  with  respect  to  Z  and 
rescaling  Z  and  Z  *  by  (2)  Substituting  the  expressions 

G(Z*)  =  f  F(Z)</Z,  =F(Z*)  (A3) 

Jo  dZ* 

into  the  inhomogeneous  integral  Eq.  (Al)  and  integrating 
by  parts  gives  the  equivalent  form: 


G(Z*)sgn(Z-Z*)erfc(lZ-Z*l)rfZ*  =  1. 

(A4) 


Differentiation  with  respect  to  Z  gives 

G(Z)  =  n  f*  G(Z*)exp[  -  (Z-Z*)-]dZ*. 

Jo 

(A5) 

If  we  now  introduce  the  change  of  variables 
Z  =  p(l)-''-,  Z*  =  t(2)-''\ 

G[p(l)  ''■]=Ag(p)  (A6) 

into  the  integral  Eq.  (A5),  then  we  obtain  Hunt  and 
Stewartson’s  integral  Eq.  (A2).  Here  ,4  is  an  arbitrary  multi¬ 
plicative  constant  since  their  integral  equation  is  homogen¬ 
eous.  The  integral  Eq.  ( Al )  has  a  unique  solution,  but  one 
must  recognize  that  F(Z)  has  a  singularity  at  Z  =  0  and 
must  adjust  the  strength  of  the  singularit;^ /„  to  obtain  a 
numerical  solution  which  gives  a  bounded  Fat  Z  =  0.  The 
integral  Eq.  ( Al )  does  not  have  a  unique  solution.  Hunt  and 
Stewartson'*'  made  the  solution  unique  with  the  condition 


dg(p) 

dp 


- 1  a.%  p- 


00, 


(A7) 


which  they  derived  by  matching  the  adjacent  core  solution. 
No  such  condition  is  needed  here,  and  F—  —  2  automatical¬ 
ly  as  Z—  00,  which  is  equivalent  to  condition  (A7). 
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